AD-A236  504 


AN  EXPERT  SYSTEM  FOR  THE  IDENTIFICATION  OF 
SILICON  IMPURITIES  USING  FTIR  SPECTROSCOPY 


Ronald  E.  Perrin 
Steven  B.  Fairchild 

Electronic  and  Optical  Materials  Branch 
Electromagnetic  Materials  &  Survivability  Division 


Final  Report  for  Period  Jan  89  to  Jun  90 


APRIL  1991 


Approved  for  Public  Release;  Distribution  Unlimited. 


MATERIALS  DIRECTORATE 

WRIGHT  LABORATORY 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT  PATTERSON  AIR  FORCE  BASE,  OHIO  45433- 


6533 


91  n  o 


NOTICE 


When  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  in  connection  with  a  definitely  Government-related  procurement, 
the  United  States  Government  incurs  no  responsibility  or  any  obligation  whatsoever. 
The  fact  that  the  government  may  have  formulated  or  in  any  way  supplied  the  said 
drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  implication,  or 
otherwise  in  any  manner  construed,  as  licensing  the  holder,  or  any  other  person  or 
corporation;  or  as  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell  any 
patented  invention  that  may  in  any  way  be  related  thereto. 


This  report  is  releasable  to  the  National  Technical  Information  Service  (NTIS). 
At  NTIS,  it  will  be  available  to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


RONALD  E.  PERRIN 
Elpbfronic  &  Optical  Materials  Branch 


STEVEN  B.  FAIRCHILD 
Hardened  Materials  Branch 


2JLnM 


ROBERT  L.  DENISON,  Chief 
Electronic  &  Optical  Materials  Branch 
Electromagnetic  Mat’ls  &  Surv.  Division 


FOR  THE  COMMANDER 


WILLIAM  R.  WOODY,  Gfiief 
Electromagnetic  Mat’ls  &  Surv.  Division 
Materials  Directorate 


If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list,  or  if 
the  addressee  is  no  longer  employed  by  your  organization  please  notify  WL/MLPO. 
WPAFB,  OH  45433-6533  to  help  us  maintain  a  current  mailing  list. 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  security 
considerations,  contractual  obligations,  or  notice  on  a  specific  document. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Id  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

lb  RESTRICTIVE  MARKINGS 

7a  SECURITY  CLASSIFICATION  AUTHORITY 

3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  Public  Release; 
distribution  unlimited. 

7b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

WL-TR-91  -4012 

5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

6a  NAME  OF  PERFORMING  ORGANIZATION 

Materials  Directorate 


li'J  il*  111  •  •!•]  r  U»]  ts 


6c  ADDRESS  (City.  State,  and  ZIP  Code) 

Wright-Patterson  AFB,  OH  45433-6533 


7a  NAME  OF  MONITORING  ORGANIZATION 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 


8a  NAME  OF  FUNDING  t  SPONSORING 
ORGANIZATION 


8c  ADDRESS  (C/fy,  State,  and  ZIP  Code) 


8b  OFFICE  SYMBOL  1  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable)  f 


10  SOURCE  OF  FUNDING  NUMBERS 


NO 


1 1  TITLE  (Include  Security  Classification) 

An  Expert  System  for  the  Identification  of  Silicon  Impurities  Using  FTIR  Spectroscopy 


17  PERSONAL  AUTHOR(S) 

Ronald  E.  Perrin  and  Steven  B.  Fairchild 


Ha  TYPE  OF  REPORT 

Final 


ESSM 

§r 

SSIBMPPiHi 

14  DATE  OF  REPORT  (Yean  Month,  Day) 

91  Apr  30 

1 7 _ COSATl  CODES _ I  18  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

FIELD  |  GROUP  |  SUBGROUP  | 

Expert-System,  electronic  materials,  impurities, 

Fourier  Transform  Infra-Red  SoectroscoDv 


19  ABSTRACT  (Continue  on  reverie  if  necessary  and  identify  by  block  number) 

An  expert-system  for  the  analysis  of  the  infra-red  spectra  of  silicon  based  electronic  materials  has 
been  developed.  The  Fourier  Transform  Infra-Red  Spectroscopy  system  produces  characteristic 
spectral  lines  for  impurities  present  in  silicon  materials  under  analysis.  A  MicroVax  workstation  with 
a  commercially  available  expert-system  shell  is  attached  to  the  spectrometer.  By  examining  these 
spectral  line  patterns  the  expert-system  is  able  to  deduce  the  identities  and  concentrations  of  the 
unknown  impurities.  The  use  of  this  program  greatly  reduces  the  analysis  time  required  for  each 
sample  and  produces  more  consistent  results. 


70  DISTRIBUTION  AVAILABILITY  OE  ABSTRACT  7 '  ABSTRACT  SECURITY  CLASSIFICATION 

(3 IJNC1  ASSIFIED'UNUM'UD  □  SAME  AS  RPT  G  DTK  uS<  RS  UHClaSSified 


77,1  NAMf  Of  RESPONSIBLE  INDIVIDUAI  77b  : E  l  E PHONE  (Include  Area  Code)  |  77c  OFFICE  SYMBOL 

rrin  or  Steve  Fairchild _  I  513/255-4474 


DDForrr  1473,  JUN  86  iv«*k.>uS  editions  .UP  obsolete  SECURITY  Cl  ASSIFICAT'ON  OF  Th! 


SECURITY  CiASS.FICAT'pN  OF  Th»S  PAGE 

UNCLASSIFIED 


TABLE  OF  CONTENTS 


SECTION  PAGE 

1 .0  INTRODUCTION  1 

1.1  BACKGROUND  2 

1.2  PROBLEM  STATEMENT  4 

1.3  PROJECT  OBJECTIVES  5 

2.0  TUTORIAL  9 

2.1  THE  REASONING  DIMENSION  1 0 

2.2  THE  REPRESENTATION  DIMENSION  1 0 

2.3  EXAMPLE  11 

3 . 0  KNOWLEDGE  BASE  DESIGN  1 7 

3.1  THE  OBJECT  NETWORK  1 9 

3.2  THE  RULE  NETWORK  23 

3.3  PROGRAM  FLOW  34 

4.0  SAMPLE  RUN  49 

5 . 0  TESTING  THE  EXPERT  SYSTEM  51 

5.1  DISCUSSION  52 

6.0  FUTURE  PLANS  56 

APPENDIX  A  FORTRAN  SUBROUTINES  57 

APPENDIX  B  LIST  OF  RULES  76 

iii 


REFERENCES 


116 


FIGURE 

1 

2 

3 

4 

5 

6 

7 

8 
9 

1  0 

1  1 

1  2 

1  3 


LIST  OF  ILLUSTRATIONS 


PAGE 


The  symbolic  representation  of  a  class  with  its  11 

related  objects  and  properties. 

A  silicon  spectrum  showing  evidence  of  the  1  3 

presence  of  the  impurity  indium. 

Symbolic  representation  of  the  class  14 

INDIUM_LINES. 

The  class  LINES.  19 

The  class  ALJJNES.  20 

The  class  LINES,  its  properties,  subclasses,  22 

and  objects. 

The  class  POSSIBLE_PEAKS.  24 

The  class  POSSIBLE_AL_PEAKS.  28 

The  rule  START_THINGS_OFF  and  the  branch  of  36 

the  rule  network  triggered  by  the  hypothesis 
PEAKS_LOADED. 

The  PEAKS_EXAMINED  branch  of  the  rule  37 

network. 

The  PEAKS_ANALYZED  branch  of  the  rule  38 

network. 

Forward  chaining  from  the  hypothesis  42 

ANALYZED. 

The  aluminum  rules.  44 


1  4  The  boron  rules. 


45 


LIST  OF  ILLUSTRATIONS 

(continued) 


FIGURE  PAGE 

1 5  The  gallium  rules.  46 

1  6  The  indium  rules.  47 

1  7  The  carbon  and  oxygen  rules.  48 

1  8  Possible  peak  boundaries  in  a  noisy  spectrum.  55 


Accession  Fop 

Ttis'craai 

DTIC  TAB 
Unannounced 
Ju ji! float lon_ 


By - - 

Distribution/ 

Availability  Codes 
j A  vail  and/or 
iDlet  j  Special 


LIST  OF  TABLES 


TABLE  PAGE 

1  Results  of  the  A!  feasibility  survey.  3 

2  All  of  the  peaks  that  can  occur  in  a  silicon  7 

spectrum  due  to  the  presence  of  the  impurities. 

3  The  report  generated  by  the  expert  system.  50 

4  Results  of  the  gallium  concentration  study  test.  53 


1.0 


INTRODUCTION 


The  Wright  Laboratory,  Electromagnetic  Materials  Division 
(WL/MLP)  conducts  research  to  develop  electronic  materials  of 
future  interest  to  the  Air  Force.  These  materials  can  include  high 
purity  silicon  and  gallium  arsenide.  These  materials  are  used  as 
infra-red  detectors  in  space  based  applications  as  well  as  being  used 
as  a  basis  material  for  high-speed  integrated  circuits  in  the  next 
generation  of  computers.  The  impurities  and  defects  in  these 
materials  determine  the  operational  limitations  of  the  devices 
constructed  from  them.  One  of  the  optical  characterization  techniques 
used  to  detect  these  defects  is  low  temperature  Fourier  Transform 
Infra-Red  spectroscopy  (FTIR).  The  material  being  examined  is  cooled 
to  a  cryogenic  temperature  and  then  excited  by  an  infra-red  light 
source.  The  absorption  and  transmittance  of  infra-red  radiation  by 
these  impurities  and  defects  produces  a  unique  spectrum.  By 
examining  the  spectrum,  an  experienced  spectroscopist  can  determine 
the  presence  of  each  impurity  as  well  as  its  concentration  in  the 
material1.  The  locations  of  peaks,  the  relative  peak  heights,  and  the 
shape  of  the  background  in  the  resulting  spectral  plot  are  used  to 
characterize  the  material's  properties.  The  presence  of  unwanted 
impurities  in  these  materials  degrades  the  performance  of  the  end 
application  devices. 
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1.1  BACKGROUND 


Modern  research  and  development  laboratories  can 
benefit  from  applying  new  automation  techniques.  The  application  of 
these  new  technologies  can  result  in  more  accurate,  more  consistent, 
and  faster  data  acquisition,  thus  maximizing  personnel  resources. 
WL/MLP  decided  that  an  Expert  System  application  to  Fourier 
Transform  Infra-Red  Spectroscopy  (FTIR)  would  be  an  optimum 
project  to  pursue  in  terms  of  Al  applications  in  the  modern 
laboratory.  In  1986  and  1987,  two  separate  groups  of  artificial 
intelligence  (Al)  experts  from  universities  and  industry2  surveyed 
the  research  being  performed  at  WL/MLP).  All  current  and  future 
areas  of  research  being  performed  within  the  WUD  48  and  WUD  50  in- 
house  laboratories  were  examined.  The  groups  were  asked  to  identify 
laboratory  systems  which  could  benefit  from  the  application  of 
artificial  intelligence  technology  and  recommend  useful  hardware 
and  software  products.  After  interviewing  the  researchers  and  being 
briefed  on  the  details  of  each  laboratory's  data  acquisition  and 
analysis  equipment,  the  Al  working  groups  generated  final  reports. 
These  reports  outlined  which  research  areas  would  benefit  from 
applying  Al  technology  and  produce  the  greatest  pay-off  to  the  Air 
Force.  This  pay-off  would  be  realized  in  terms  of  higher  productivity 
and  more  accurate  data  analysis  if  automated  using  Al  technology. 
The  results  of  one  survey  are  listed  in  Table  1.  Note  that  the  areas 
referenced  as  FTIR-2  and  FTIR-3  show  the  highest  ranking  in  terms 
of  feasibility  and  pay-off.  We  proposed  that  an  expert  system  that 
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Table  1.  Results  of  Al  feasibility  survey. 


OVERALL  CONCLUSIONS  AND  RECOMMENDATIONS 
The  rankings  given  in  the  table  for  the  various  Al  application  opportunities  are  based  on 


technical  considerations  alone.  (  1  =  highest  potential  ) 

Opportunity  Number  and  Name 

Recommendation 

Rank 

MBE-1 

Utilize  Past  Experience 

pursue/guarded/ 
build  data-base  first 

3 

MBE-2 

Interp.  RHEED  Oscill. 

Do  not  pursue 
re-examine  later 

MBE-3 

Interp.  Optical  Probes 

Do  not  pursue 
re-examine  later 

MBE-4 

Qual.  Simul.  of  Growth 

Pursue/guarded/ 
build  data-base  first 

4 

PVD-1 

Correction  Analysis 

Do  not  pursue 

PVD-2 

Detection  Analysis 

Do  not  pursue 

PVD-3 

Intell.  Front-End  to 

Chem.  Abstracts  Data-Base 

Pursue  with  constraints 

2 

FTIR-1 

Initial  Set  of  Experiments 

Do  not  pursue 

FTIR-2 

ID  Known  Centers  and  Locate 

Unknown  Centers 

Pursue 

1 

FTIR-3 

Categorize  and  Assess  Significant 

Unknown  Centers 

Pursue 

2 

FTIR-4 

Correlate  Unknown  Centers 

Pursue/guarded 
build  data-base  first 

3 

FTIR-5 

Formulate  Strategy  for  Further 

Experiments 

Pursue/guarded/ 
build  data-base  first 

4 

3 


could  address  the  areas  described  in  Table  1  as  FTIR-2  and  3  be 
designed.  This  project  would  not  only  benefit  the  FTIR  laboratory  but 
lessons  learned  would  be  extremely  valuable  for  future  Al  projects 
within  WL/MLP. 

1.2  PROBLEM  STATEMENT 

The  task  of  analyzing  FTIR  spectra  is  complex  and  tedious. 
Simply  matching  peaks  with  known  impurity  wavenumber  locations 
tells  the  researcher  very  little.  The  most  common  impurities, 
aluminum,  boron,  gallium,  indium,  carbon  and  oxygen,  all  occur  in  the 
same  area  of  the  spectrum  (0-2000  wavenumbers).  The  presence  of 
a  singularity,  or  peak,  in  the  spectrum  at  a  certain  wavenumber 
could  be  the  result  of  dust  present  in  the  sample  chamber  or  it  could 
indicate  the  presence  of  an  impurity  or  sample  defect.  The  thickness 
of  the  sample  as  well  as  the  method  of  mounting  it  also  affects  the 
shape  of  the  spectrum.  Sixty-hertz  pickup  and  its  harmonics  show  up 
in  spectra  as  well  as  noise  spikes  that  occasionally  occur  due  to 
power  line  transients.  The  researcher  must  carefully  look  for  peak 
patterns  knowing  from  experience  that  many  false  conclusions  could 
be  reached  if  certain  rules  of  thumb  are  not  followed.  It  is  this  well 
defined  set  of  rules  and  heuristics  which,  if  incorporated  into  a 
knowledge  base,  could  allow  a  less  experienced  researcher  or  aide  to 
accurately  characterize  samples.  This  type  of  knowledge  base  could 
be  useful  for  basic  research  as  well  as  for  routine  material  quality 
control  testing. 
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Since  the  area  of  a  selected  peak  is  used  to  calculate  its 
impurity  concentration,  the  expert  system  must  be  capable  of 
defining  the  baseline  of  a  peak.  While  this  task  may  seem  simple,  it 
must  be  remembered  that  peaks  can  be  very  weak  if  the 
concentrations  are  small.  Peaks  can  also  occur  in  very  r^isy  areas 
of  the  spectrum  or  may  be  present  where  the  spectral  background  is 
increasing  or  decreasing.  Choosing  the  baseline  of  a  peak  in  a  noisy 
area  is  a  subjective  choice  based  on  the  past  experience  of  the 
researcher.  Often  two  expert  spectroscopists  will  calculate  peak 
areas  differing  by  30%  or  more  when  examining  the  same  peak  in  a 
noisy  background.  The  process  of  calculating  peak  areas  accurately 
can  be  extremely  time  consuming.  Researchers  have,  in  the  past, 
resorted  to  cutting  the  peak  of  interest  out  of  a  spectral  plot  and 
weighing  it  on  a  set  of  sensitive  scales.  This  weight  was  compared 
to  the  weight  of  a  unit  area  from  the  same  type  of  graph  paper. 
Obviously  an  easier,  more  consistent  method  is  needed.  The 
application  of  the  proposed  expert-system  would  generate  more 
consistent  analyses  and  make  a  tedious  and  exacting  task  quick  and 
routine. 

1.3  PROJECT  OBJECTIVES 

One  of  the  most  difficult  tasks  in  designing  this  project  was 
deciding  what  the  expert  system  should  actually  be  able  to  do.  There 
is  always  a  trade-off  between  what  is  needed  by  the  end  user  and 
what  is  feasible  to  design  given  the  constraints  of  of  the  expert- 
system  shell,  the  skills  of  the  designer,  and  a  reasonable  completion 
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date.  There  is  a  saying,  "Anything  can  be  made  to  work  if  you  tinker 
with  it  long  enough."  We  wanted  to  limit  the  scope  of  this  project  to 
be  realistic  given  our  expert  system  design  experience,  time 
constraints,  and  our  other  job  responsibilities. 

It  was  decided  that  the  expert  system,  given  an  FTIR  absorption 
coefficient  spectrum,  should  be  able  to  identify  all  impurities 
present,  their  concentrations  in  the  material,  and  all  the  associated 
impurity  lines  present.  If  the  spectrum  is  a  photoconductivity 
spectrum,  it  should  be  able  to  identify  all  impurities  present  and 
their  associated  line  numbers. 

To  accomplish  this,  a  knowledge  base  had  to  be  designed  that 
incorporated  all  the  necessary  information  about  how  to  identify 

these  impurities  in  a  silicon  spectrum.  The  theoretical  background 

for  this  knowledge  base  came  from  two  sources.  The  first  source  is 
the  knowledge  obtained  from  extensive  question  and  answer  sessions 
with  experienced  spectroscopists  in  this  particular  area  of  materials 
analysis.  This  part  is  typically  refered  to  as  'knowledge  engineering.' 
The  second  source  of  information  is  technical  papers  and  reports 
directed  in  this  area  of  research.  Together  these  sources  provide  a 
realistic  starting  point  for  constructing  a  useful  expert  system. 
Table  2  is  taken  from  such  a  paper2  and  lists  all  the  lines  (i.e. ,  peaks) 

by  wavenumber  that  can  occur  in  a  silicon  spectrum  owing  to  the 

impurities  of  interest. 

There  should  be  minimal  user  interaction.  The  user  could  be 
asked  if  they  agree  with  interim  results;  if  not,  the  user  has  the 
authority  to  modify  these,  but  the  final  conclusions  are  determined  by 
the  expert  system.  A  graph  of  the  spectrum  being  examined  should  be 
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Table  2.  All  of  the  peaks  that  can  occur  in  a  silicon  spectrum 
due  to  the  presence  of  the  impurities. 


Line 

BORON 

ALUMINUM 

r  *  W 

INDIUM 

THALLIUM 

1 

244.96 

442.73 

469.75 

1145.41 

1876.1 

2 

278.19 

471.93 

— 

1175.82 

1906.7 

3 

309.32 

— 

541.45 

1208.00 

1937.4 

4 

319.32 

516.52 

548.23 

1216.31 

1947.9 

4a 

319.97 

524.11 

550.57 

1218.57 

4b 

321.9 

525.56 

552.02 

1219.46 

5 

334.49 

535.16 

564.2 

1232.2 

6a 

338.03 

— 

— 

6b 

339.24 

— 

567.82 

— 

6 

340.05 

538.71 

568.71 

1236.56 

1968.0 

7 

342.07 

541.21 

570.97 

1239.05 

8 

344.73 

543.87 

573.72 

1241.96 

8a 

346.18 

545.56 

575.49 

1243.33 

9a 

348.12 

546.94 

576.94 

— 

9 

348.76 

547.74 

577.67 

1245.83 

10a 

350.78 

549.76 

579.72 

1247.93 

10b 

351.99 

550.24 

— 

10 

352.55 

557.53 

581.30 

1249.62 

10c 

353.12 

552.74 

581.78 

lOd 

354.49 

553.47 

583.23 

1250.83 

11a 

— 

554.92 

584.28 

— 

11b 

356.5 

555.49 

584.93 

1253.41 

11 

357.07 

556.21 

585.73 

1254.3 

Ei 

358.04 

556.78 

586.62 

1254.79 

1985.8 

2p' 

668.57 

867.31 

897.23 

1565.4 

2285.0 

3p‘ 

692.85 

891.75 

921.51 

1589.84 

2309.2 

4p' 

701.32 

900.22 

929.98 

1598.47 

5p' 

— 

— 

934.01 

Ei* 

712.93 

911.75 

941.5 

1609.76 

XI 

184 

338 

347 

799 

X2 

220 

372 

381 

831 

X3 

250 

404 

— 

863 

X4 

261 

414 

423 

873 

Ei(x) 

299 

454 

460 

909 

displayed  with  all  peaks  marked  by  the  expert  system.  If  any  peaks  of 
interest  have  been  missed,  the  user  should  be  able  to  add  these  with 
the  mouse  and  cursor.  The  user  should  not  be  allowed  to  delete  peaks 
that  the  system  has  marked;  it  is  up  to  the  expert  system  to  conclude 
whether  marked  peaks  are  significant.  All  impurities  identified  must 
be  present  in  the  patterns  that  experienced  spectroscopists  recognize 
as  being  necessary  to  insure  a  high  degree  of  confidence.  Noise  and 
spurious  peaks  present  in  the  spectrum  should  be  ignored. 
Indeterminate  peaks  should  be  noted  and  classified  as  unknown. 
Finally,  a  report  containing  this  information  should  be  prepared  in  a 
concise  format  and  be  displayed  and  printed  upon  conclusion  of  the 
analysis. 
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2.0  TUTORIAL 


The  purpose  of  an  Al  knowledge  based  system  is  to  transform 
the  line  of  thought  of  a  human  expert  into  a  set  of  formalized 
symbolic  operations  so  as  to  reproduce  his  or  her  expertise.3 
Knowledge  design  is  the  process  of  transferring  this  formulated 
knowledge  from  humans  to  the  machine.  Once  the  computer  has  this 
set  of  rules  formulated  from  human  logic,  it  can  perform  deductive 
or  inductive  reasoning  and  arrive  at  the  same  conclusion  as  the 
human  expert  when  given  the  appropriate  data. 

The  methodology  involved  in  knowledge  engineering  involves 
two  closely  related  critical  dimensions:  the  reasoning  dimension 
and  the  representation  dimension .3  The  reasoning  dimension  is 
composed  of  the  rules  that  constitute  the  logic  the  expert  system 
uses  to  perform  reasoning.  The  representation  dimension  is 
composed  of  objects  and  classes,  which  can  be  thought  of  as 
symbolic  representations  of  the  data  that  the  rules  reason  upon. 
Therefore,  objects  and  classes  are  complimentary  to  rules. 

The  reasoning  dimension  and  the  representation  dimension  are 
closely  interrelated  to  provide  an  efficient  network.  They  allow 
ease  in  knowledge  integration,  which  is  the  reduction  of  thought 
mechanisms  to  more  understandable  concepts.  Once  knowledge 
integration  has  been  accomplished,  powerful  symbolic  computations 
can  be  performed  which  allow  automation  of  the  problem  solving 
process.3 
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2.1  THE  REASONING  DIMENSION 


To  represent  reasoning,  an  expert  system  uses  rules.  Rules  are 
symbolic  structures  which  express  deductive  or  evocative 
progressions  in  a  reasoning  path.  By  deductive  we  mean  that  a  rule 
can  be  used  to  verify  conditions  in  another  rule.  By  evocative,  we 
mean  that  a  rule  can  trigger  the  activation  or  the  evaluation  of  other 
rules.3 

A  rule  is  a  chunk  of  knowledge  representing  a  situation  and  its 
immediate  consequences.  The  format  of  a  rule  is  a  symbolic 
structure  of  the  type 

if  ...  then  ...  and  do... 

where  if  is  followed  by  a  set  of  conditions,  then  by  a  hypothesis  or 
goal  which  becomes  TRUE  when  the  conditions  are  met,  and  do  by  a 
set  of  actions  to  be  undertaken  as  a  result  of  a  positive  evaluation 
of  the  conditions  of  the  rule.3 

2.2  THE  REPRESENTATION  DIMENSION 

The  representation  dimension  is  composed  of  a  set  of 
interrelated  objects  which  represent  the  data  or  events  that  the 
rules  act  upon.  Since  the  rules  reason  on  a  representation  of  the 
problem,  the  reasoning  and  representation  dimension  are  inherently 
related.3 

The  representation  of  knowledge  is  accomplished  with  three 
symbolic  tools:  objects,  properties,  and  classes.  Various 


combinations  of  these  three  components  provide  an  efficient  method 
of  representating  the  problem  upon  which  the  rules  act.3 

An  object  is  an  elementary  unit  of  description.  Everything  is 
an  object.  A  property  is  a  characteristic  of  an  object.  A  class  is  a 
collection  of  objects  that  share  properties.  As  illustrated  in  Figure 
1,  objects  are  illustrated  by  triangles,  properties  by  squares,  and 
classes  by  circles.3 


Figure  t.  The  symbolic  representation  of  a  class  with  its 
related  objects  and  properties.3 

2.3  EXAMPLE 


This  relationship  between  classes,  objects,  and  properties  and 
their  connection  to  rules  can  be  seen  in  the  following  example, 
which  introduces  the  concept  of  knowledge  base  design  to 
spectroscopic  analysis. 


A  spectroscopist  can  examine  the  optical  spectrum  of  a 
material  to  determine  its  composition.  For  example,  pure  silicon 
can  always  be  identified  by  its  specific  characteristic  spectrum.  If 
a  certain  spectrum  of  silicon  contains  peaks  that  would  not  show  up 
in  a  spectrum  of  pure  silicon,  then  there  are  impurities  present  in 
the  silicon  sample  from  which  the  spectrum  was  generated.  From 
previous  experimental  knowledge  we  know  that  various  impurities 
identify  themselves  by  the  specific  wavenumbers  at  which  these 
peaks  appear. 

For  example,  if  the  impurity  indium  is  present  in  a  silicon 
sample,  it  can  manifest  itself  by  causing  as  many  as  20  additional 
peaks  to  appear  in  the  spectrum,  three  of  which  can  be  used  to 
confirm  the  presence  of  indium.  These  peaks  will  be  referred  to  as 
lndium_Lines.  The  confirming  peaks  are  lndium_Line_1 , 
lndium_Line_2,  and  lndium_Line_4.  All  three  of  these  lines  must  be 
present  for  indium  to  exist.  Figure  2  shows  an  example  of  a  peak 
pattern  in  a  silicon  spectrum  that  identifies  the  presence  of  indium. 

Figure  3  illustrates  the  object  structure.  lndium_Lines  is  a 
class,  and  lndium_Line_1 ,  2  &  4  are  the  objects.  As  previously 
stated,  a  class  is  a  group  of  objects  that  have  the  same  properties. 
Since  each  object  lndium_Line_n  wili  be  identified  by  a  specific 
wavenumber,  each  object  will  be  given  the  property  WN,  for 
wavenumber.  Each  object  will  also  be  given  the  property  EXIST, 
which  will  be  used  to  identify  whether  or  not  that  particular  peak  is 
present. 


1e  FTIR  Absorbance  Spectrum 


Figure  2.  A  silicon  spectrum  showing  evidence  of  the  presence  of  i 


Indium-Lines 


Figure  3.  Symbolic  representation  of  the  class  lndium_Lines. 

Rules  must  operate  on  data,  and  data  in  the  object  network  is 
represented  by  slots.3  A  slot  is  the  value  that  is  held  in  the  property 
of  an  object: 

object. property  =  slot 

For  this  object  network,  a  slot  would  be  lndium_Line_1  .EXIST  and 
this  slot  (datum)  would  have  a  boolean  value.  Another  slot  would  be 
lndium_line_4.WN,  and  it  would  have  a  numeric  value.  For  example, 
the  wavenumber  for  the  third  indium  line  is  1216,  so  the  value  held 
in  the  slot  lndium_line_4.WN  is  1216. 

Rules  operate  on  the  objects,  which  are  a  symbolic 
representation  of  the  problem.3  A  rule  can  be  designed  that  will 
confirm  whether  or  not  indium  is  present  in  the  silicon  spectrum. 
This  rule  attempts  to  confirm  the  presence  of  the  confirmation 
pattern. 


if  lndium_Line_1  .EXIST  is  true 

and  lndium_Line_2. EXIST  is  true 

and  lndium_Line_4. EXIST  is  true 

then  lndium_Preser>t  is  confirmed 

In  this  rule,  there  are  three  conditions  which  lead  to  the 

evaluation  of  a  hypothesis.  The  conditions  are  the  tests  for  the 

presence  of  the  three  indium  lines,  and  the  hypothesis  is  the  attempt 
to  confirm  the  presence  of  indium.  All  three  conditions  must  be  true 

in  order  for  the  hypothesis  to  be  confirmed,  i.e.,  all  three  indium 

lines  must  be  present  in  order  for  the  presence  of  indium  to  be 
confirmed. 

This  rule  operates  on  the  slots  that  were  set  up  in  the  object 
network.  It  will  evaluate  the  slots  lndiumjine_n.  EXIST.  (Assume 
that  a  previous  rule  has  determined  the  values  for 
lndium_line_n.WN.)  The  values  of  these  slots  will  determine 
whether  or  not  the  hypothesis  is  confirmed.  The  name  of  the 
hypothesis  is  lndium_  Present,  which  has  a  Boolean  value.  If  all  the 
data  values  (i.e.,  lndium_line_/7. EXIST)  that  occur  in  the  conditions 
evaluate  to  true,  then  the  hypothesis  lndium_Present  will  evaluate 
to  true.  If  any  of  the  datum  evaluate  to  false,  the  hypothesis 
lndium_  Present  will  evaluate  to  false. 

Therefore,  if  the  pattern  of  the  three  indium  lines  is  found, 
this  rule  will  confirm  that  indium  is  indeed  present.  If  any  one  of 
the  three  lines  does  not  exist,  then  the  rule  confirms  that  indium  is 
not  present,  since  the  whole  pattern  must  be  present  to  confirm  the 
presence  of  indium. 


This  example  has  introduced  the  concept  of  knowledge  base 
design.  The  most  important  part  of  knowledge  base  design  is  the 
inherent  relationship  between  rules  and  the  object  network.  As 
illustrated,  the  object  network  is  composed  of  classes,  objects,  and 
properties,  which  are  a  symbolic  representation  of  the  problem. 
Rules  can  be  related  to  the  object  network  through  their  ability  to 
evaluate  data.  Data  are  represented  by  a  slot,  which  is  the  value 
stored  in  a  property  attached  to  an  object. 


3.0  KNOWLEDGE  BASE  DESIGN 


As  previously  stated,  the  critical  task  of  this  project  is  to 
design  the  knowledge  base  so  that  the  expert  system  will  have  the 
ability  to  identify  all  impurities  of  interest  that  can  be  present  in  a 
silicon  sample.  These  impurities  are  aluminum,  boron,  galium, 
indium,  carbon,  and  oxygen.  The  presence  of  these  impurities  can  be 
confirmed  by  identifying  the  variations  they  cause  in  a  silicon 
absorption  spectrum.  This  confirmation  process  involved  three 
critical  steps:  (1)  defining  the  wavenumber  range  where  the  peaks 
occur  for  each  impurity,  (2)  identifying  the  patterns  that  confirm 
their  existence,  (3)  and  classifying  all  additional  lines  associated 
with  that  impurity  which  are  not  contained  within  the  identifying 
pattern. 

For  the  first  step,  the  wavenumber  ranges  of  the  various 
impurities  must  be  identified.  From  previous  experimental 
knowledge,  it  is  known  that  each  possible  impurity  causes  peaks  to 
occur  in  a  well  defined  region.  For  example,  if  aluminum  is  present 
in  a  silicon  sample,  peaks  will  appear  in  the  spectrum  in  a  range 
from  337  to  912  wavenumbers.  As  many  as  20  different  peaks  can 
show  up  in  this  range.  The  problem  lies  in  determining  which  peaks 
present  in  this  range  are  actually  due  to  the  presence  of  aluminum 
and  which  ones  are  caused  by  extraneous  noise. 

This  leads  to  the  second  step,  identifying  the  patterns  of  each 
impurity  to  confirm  its  presence.  Each  of  the  20  aluminum  peaks 
that  may  appear  can  be  distinctly  defined  by  a  particular 
wavenumber,  which  we  will  call  aluminum_line_n.  If  aluminum  is 


present,  it  is  known  from  previous  experimental  knowledge  that  the 
pattern  of  aluminumjines  1,  2,  and  8  must  appear  in  the  spectrum. 
If  strong  identifiable  peaks  appear  around  442  wavenumbers 
(aluminumjine_1),  471  wavenumbers  (a!uminum_line_2),  and  543 
wavenumbers  (aluminum_iine_8),  then  you  can  positively  conclude 
that  aluminum  is  present. 

After  this  pattern  matching  has  accomplished  the 
identification  of  the  impurity,  the  next  step  is  to  identify  all 
additional  peaks  that  are  associated  with  that  impurity.  For 
aluminum,  there  may  be  as  many  as  17  peaks,  since  only  3  out  of  a 
possible  20  were  needed  to  identify  the  pattern.  These  peaks  must 
be  identified  by  their  wavenumbers  and  classified  as 
aluminum_line_n.  The  important  point  to  remember  is  that  there 
may  be  peaks  that  show  up  at  the  aluminum_line_n  wavenumbers, 
that  are  not  due  to  the  presence  of  aluminum.  If  the  pattern  of 
aluminumjines  1,  2,  and  8  has  not  been  previously  confirmed,  the 
presence  of  these  peaks  must  be  caused  by  some  extraneous  noise. 
Any  peaks  that  show  up  at  the  aluminumJine_/7  wavenumbers  in  the 
absence  of  the  confirmation  of  the  aluminum  pattern  would  be 
considered  unidentified. 

In  summary,  the  goal  of  this  effort  is  to  accomplish  the 
impurity  analysis  for  the  six  possible  impurities  of  interest  that 
can  occur  in  silicon.  These  are  aluminum,  boron,  gallium,  indium, 
carbon  and  oxygen.  The  identification  process  is  the  same  for  each 
impurity.  First,  the  wavenumber  range  where  they  will  cause  peaks 
to  occur  must  be  defined.  Next,  the  patterns  of  peaks  that  confirm 
their  existence  must  be  identified.  And  finally,  all  additional  peaks 


that  are  caused  by  the  impurity  which  are  not  contained  in  the 
identifying  pattern  must  be  classified. 

3.1  THE  OBJECT  NETWORK 


As  stated  in  the  tutorial  section,  the  object  network 
represents  the  representation  dimension.  The  representation  is 
composed  of  classes,  objects,  and  properties  which  can  be  thought 
of  as  symbolic  representations  of  the  data  that  the  rules  reason 
upon.  The  object  network  designed  for  this  situation  must  therefore 
use  classes,  objects,  and  properties  to  represent  the  task  of  using 
spectral  analysis  to  identify  impurities  in  silicon. 

The  peaks  that  appear  in  a  silicon  spectrum  that  are  due  to  the 
presence  of  impurities  are  classified  as  lines.  Therefore  a  class 
will  be  used  called  LINES.  It  will  have  six  subclasses  for  the  six 
impurities  that  may  be  .  "esent,  aluminum,  boron,  gallium,  indium, 
carbon,  and  oxygen  (Figure  4).  These  subclasses  will  be  called 
AL_LINES,  BOJJNES,  GAJJNES,  INJJNES,  CAJJNES,  and  OXJJNES. 
Classes  are  symbolically  represented  by  circles. 


o  LINES 


O  AL  LINES 
O  BOJJNES 

Oga_lines 
O  IN  LINES 
O  CA  LINES 
O  OX  LINES 


Figure  4.  The  class  LINES 


Each  subclass  will  contain  objects  to  represent  all  of  the  lines 
whose  existence  may  be  due  to  the  presence  of  that  impurity.  For 
example,  the  subclass  AL_LINES  (Figure  5)  will  contain  the  objects 
AL_LINE_1,  AL_LINE_2,  etc.  Each  one  of  these  objects  has  two 
properties,  EXIST,  and  WN  (wavenumber).  EXIST  will  have  a  Boolean 
value  and  WN  will  have  a  numeric  value.  The  slot  AL_LINE__1  EXIST 
will  be  set  to  TRUE  if  that  particular  line  exists.  If  the  line  exists, 
the  slot  AL_LINE_1.WN  will  be  given  the  value  of  the  wavenumber 
where  the  peak  was  found.  Objects  are  symbolically  represented  by 
triangles  and  properties  by  squares. 
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In  the  previous  diagram,  the  properties  EXIST  and  WN  are 
shown  attached  to  the  class  ALJJNES,  rather  than  to  the  individual 
objects.  This  illustrates  the  principle  of  downward  inheritance.  If 
properties  are  initially  attached  to  a  class,  then  all  the  objects 
within  that  class  will  inherit  these  properties.  This  is  very 
efficient  for  the  situation  where  all  of  the  objects  in  a  class  are 
similar  items  with  the  same  characteristics.  The  properties  EXIST 
and  WN  don't  have  to  be  individually  attached  to  each  of  the  20 
aluminum  lines.  As  long  as  they  are  attached  to  the  class,  these 
properties  will  be  inherited  by  all  20  objects  which  represent  the 
aluminum  lines. 

Figure  6  shows  the  full  extent  of  the  class  LINES.  It  is  shown 
with  all  of  its  subclasses,  and  all  of  the  objects  that  are  contained 
within  each  subclass.  As  shown,  downward  inheritance  also  occurs 
from  class  to  subclass.  Initially,  the  properties  EXIST  and  WN  are 
attached  to  the  class  LINES.  They  are  then  inherited  down  by  the 
subclasses  ALJJNES,  BOJJNES,  etc.,  and  from  there  they  are 
inherited  down  to  each  individual  object  of  all  of  the  subclasses. 

The  objects  illustrated  represent  all  of  the  lines  that  the 
expert  system  will  attempt  to  identify.  Aluminum,  boron,  gallium, 
and  indium  are  responsible  for  the  majority  of  the  impurity  induced 
peaks  that  can  occur  in  a  silicon  spectrum.  Carbon  and  oxygen  are 
responsible  for  causing  a  peak  to  occur  at  only  one  wavenumber  each 
if  they  are  present. 
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The  presence  of  impurities  in  a  sample  of  silicon  may  cause  as 
many  as  84  different  peaks  to  appear  in  its  spectrum.  The  class 
LINES  contains  six  subclasses  which  contain  all  of  the  lines  that  are 
potentially  caused  by  these  impurities. 

3.2  THE  RULE  NETWORK 


The  rule  network  represents  the  reasoning  dimension.  The 
reasoning  dimension  is  composed  of  the  rules  that  constitute  the 
logic  the  expert  system  uses  to  perform  reasoning.  Rules  were 
designed  to  categorize  the  peaks  according  to  the  wavenumber 
ranges  of  the  impurities  of  interest,  and  also  to  identify  their  other 
lines  if  their  patterns  were  recognized. 

First,  the  expert  system  must  obtain  a  list  of  the  peaks  that 
appeared  on  the  silicon  spectrum.  After  the  spectrum  has  been 
generated  on  the  FTIR  spectrometer,  a  program  is  run  which  locates 
the  peaks  in  the  spectrum,  and  lists  them  in  a  table  by  wavenumber 
and  amplitude.  This  peak  table  has  the  following  format: 


Peak_1 

432.56 

0.654 

Peak_2 

458.78 

0.421 

Peak  3 

543.89 

0.784 

If  n  peaks  were  found  in  the  spectrum,  n  peaks  will  be  listed  in  this 
peak  table  in  ascending  order  by  wavenumber.  Peak_1  would  have 
the  lowest  wavenumber  and  Peak_n  the  highest. 
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This  peak  table  is  read  into  the  knowledge  base.  A  class  is 
created  called  POSSIBLE_PEAKS,  and  each  peak  is  an  object  of  this 
class  (Figure  7).  The  name  of  the  objects  will  be  PEAK_1  through 
PEAK_n,  and  they  will  have  the  properties  WAVENUMBER  and 
AMPLITUDE.  Using  the  previous  table  as  an  example,  the  slot 
PEAK_1  .WAVENUMBER  would  contain  the  value  432.56,  and  the  slot 
PEAK_1  .AMPLITUDE  would  contain  the  value  0.654.  With  this 
symbolic  representation  of  the  initial  peak  data,  the  rules  can  now 
begin  the  reasoning  process. 


The  first  task  of  the  rule  network  is  to  categorize  the  peaks  in 
terms  of  what  impurity  they  might  possibly  represent.  Since  the 
wavelength  range  for  aluminum  peaks  is  337  to  912  wavenumbers,  a 
rule  must  be  written  that  finds  all  of  the  peaks  in  the  class 


POSS!BLE_PEAKS  that  fall  within  this  wavelength  range.  The  name 
of  this  is  PICK_POSSIBLE_ALUMINUM  _PEAKS. 

This  rule  will  utilize  a  technique  called  pattern  matching, 
which  permits  the  creation  of  a  list  of  objects  that  satisfy  a  certain 
set  of  conditions.  First,  the  rule  will  sort  through  all  the  peaks  in 
the  class  POSSIBLE_PEAKS  and  keep  the  ones  whose  wavenumber  is 
greater  than  337.  Then  it  will  sort  through  that  list  again  and  keep 
all  of  the  peaks  whose  wavenumber  is  less  than  912.  Then  it  will 
take  this  final  list  and  use  it  to  create  a  new  class  called  POSSIBLE 
_AL_PEAKS.  This  new  class  will  contain  all  of  the  peaks  whose  peak 
location  falls  within  the  range  of  337  to  912  wavenumbers.  The  rule 
looks  as  follows: 

Name:  PICK_POSSIBLE_ALUMINUM_PEAKS 

Conditions: 

>  <|POSSIBLE_PEAKS|>.WAVENUMBER  337 

<  <|POSSIBLE_PEAKS|>.WAVENUMBER  912 

Hypothesis:  POSSIBLE_ALUMINUM 

Actions: 

CreateObject  <|POSSIBLE_PEAKS|>  |POSSIBLE_AL_PEAKS| 

As  seen  above,  rules  are  composed  of  three  parts:  conditions,  a 
hypothesis,  and  actions.  Conditions  are  conditional  statements 
which  must  be  evaluated  to  a  Boolean  true  or  false  value.  If  all  of 
the  conditions  evaluate  to  true,  the  hypothesis,  which  is  also  a 
Boolean  value,  will  evaluate  to  true,  and  the  actions  will  be 
performed.  If  even  one  of  the  conditions  evaluates  to  false,  the 
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hypothesis  will  evaluate  to  false,  and  the  actions  will  not  be 
performed.  In  the  preceding  rule,  there  are  two  conditions  and  one 
action.  If  both  of  the  conditions  evaluate  to  true,  the  hypothesis 
POSSIBLE_ALUMINUM  will  evaluate  to  true,  and  the  action  will  be 
performed.  If  the  first  condition  evaluates  to  false,  the  hypothesis 
will  automatically  evaluate  to  false  and  the  second  will  not  even  be 
evaluated. 

The  vertical  bars  enclosing  POSSIBLE_PEAKS  indicate  that  it 
is  a  class.  The  angle  brackets  outside  the  vertical  bars  indicate  that 
pattern  matching  will  occur.  The  attached  property,  WAVENUMBER, 
indicates  that  the  WAVENUMBER  slot  of  each  object  in  the  class 
POSSIBLE_PEAKS  will  be  used  for  the  comparison  that  is  determined 
by  the  conditions. 

The  first  condition  tells  the  expert  system  to  pattern  match, 
which  means  to  look  at  each  object  (PEAK_1  through  PEAK_40)  in 
the  class  POSSIBLE_PEAKS,  and  compare  the  value  of  the  slot 
PEAK_n. WAVENUMBER  to  337.  If  the  value  of  this  slot  is  greater 
than  337,  the  object  PEAK_n  is  kept.  If  it  finds  at  least  one  match, 
the  condition  will  evaluate  to  true  and  a  list  will  remain  which 
contains  all  the  objects  whose  slot  PEAK_n. WAVENUMBER  is  greater 
than  337.  If  it  doesn’t  find  a  match,  the  condition  will  evaluate  to 
false,  the  hypothesis  POSSIBLE_  ALUMINUM  will  be  set  to  false  and 
no  further  conditions  will  be  evaluated.  This  would  indicate  that 
none  of  the  peaks  in  the  class  POSSIBLE_PEAKS  has  a  peak  location 
wavenumber  that  is  greater  than  337,  so  they  can't  possibly  be 
aluminum. 
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The  second  condition  only  gets  evaluated  if  the  first  condition 
actually  creates  a  list  of  objects  whose  PEAK_n. WAVENUMBER  slot 
is  greater  than  337.  It  says  to  analyze  this  list  of  objects  and  keep 
the  ones  whose  PEAK_n.WAVENUMBER  slot  is  less  than  912.  If  it 
finds  at  least  one  match,  the  condition  will  evaluate  to  true  and  a 
list  will  remain  which  contains  all  the  objects  whose 
PEAK_n.WAVE-  NUMBER  slot  is  greater  than  337  and  less  than  912. 
If  it  doesn't  find  a  match,  this  condition  will  evaluate  to  false  which 
causes  the  hypothesis  POSSIBLE_ALUM!NUM  to  be  set  to  false.  This 
indicates  that  no  peaks  were  found  that  fall  in  the  wavenumber 
range  established  by  the  presence  of  aluminum. 

If  both  conditions  evaluate  to  true,  the  actions  will  be 
performed.  A  list  of  peaks  will  remain  which  might  possibly  be 
aluminum  since  they  fall  within  its  wavenumber  range.  Therefore,  a 
new  class  will  be  created  called  POSSIBLE_AL_PEAKS.  As  seen  in 
the  example  rule,  this  class  will  be  created  with  the  CreateObject 
command.  The  second  term  after  the  CreateObject  command, 
|POSSIBLE_AL_PEAKS|,  tells  the  expert  system  to  create  the  new 
class.  The  first  term  after  the  CreateObject  command, 
<|POSSIBLE_PEAKS|>,  represents  the  list  that  was  created  by  the 
first  two  conditions.  So,  the  CreateObject  command  takes  the  list 
of  objects  that  was  compiled  by  virtue  of  the  first  two  conditions 
and  puts  it  into  the  new  class  POSSIBLE_AL_PEAKS. 

For  example,  say  the  class  POSSIBLE_PEAKS  contains  100 
peaks  (PEAK_1  through  PEAK_100),  and  the  rule 
PICK_POSSIBLE_ALUMINUM_PEAKS  found  15  peaks  that  fall  within 
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the  aluminum  wavenumber  range.  The  class  POSSIBLE_  AL_PEAKS 
would  then  contain  15  peaks  (Figure  8). 

The  next  step  is  to  analyze  the  15  peaks  in  the  class 
POSSIBLE_AL_PEAKS  and  determine  if  the  pattern  exists  which 
confirms  the  presence  of  aluminum.  If  three  of  these  peaks  have 
wavenumbers  corresponding  to  aluminum  lines  1,  2  and  8,  then  it  can 
be  concluded  that  aluminum  is  present.  The  wavenumbers  for 
aluminum  lines  1,2,  and  8  are  approximately  442.7,  471.9,  and  443.9. 


WAYENUMBER 

AMPLITUDE 


O  POSSIBLE_AL_PEAKS 


PEAK_1 
PEAK_2 
PEAK_8 
PEA K_9 
PEAK_1 1 
PEAK_1  5 
PEAK  18 
PEAK_1 9 
PEAK_22 
PEAK  23 
PEAK  29 
PEAK  33 
PEAK  36 
PEAK_38 

PEAK_40 


Figure  8.  The  class  POSSIBLE_AL_PEAKS 


A  rule  must  be  written  that  looks  at  the  fifteen  slots 
PEAK_n. WAVENUMBER  to  determine  if  there  is  a  wavenumber  match 
with  aluminum  lines  1,  2,  and  8.  The  name  of  this  rule  will  be 
INVESTIGATE_ALUMINUM.  It  appears  on  the  following  page. 
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The  first  condition  evaluates  the  status  of  POSSIBLE_ 
ALUMINUM,  which  is  the  hypothesis  of  the  rule  ’PICK_POSSIBLE_ 
ALUMINUM_PEAKS.'  The  operator  YES  evaluates  a  Boolean  value.  This 
is  equivalent  to  asking  "Is  the  value  of  the  hypothesis  POSSIBLE_ 
ALUMINUM  true?"  If  POSSIBLE_ALUMINUM  is  true,  this  means  that 
the  previous  rule  found  some  peaks  that  might  possibly  be  aluminum 
and  stored  them  in  the  class  POSSIBLE_AL_PEAKS.  These  peaks  then 
need  to  be  analyzed  to  see  if  the  confirming  pattern  exists.  This  is 
accomplished  with  the  rest  of  the  conditions.  If  the  value  of  the 
hypothesis  POSSIBLE_ALUMINUM  is  false,  no  possible  aluminum 
peaks  were  found,  and  the  class  POSSIBLE_AL_PEAKS  was  never 

Name:  INVESTIGATE_ALUMINUM 
Conditions: 


YES 

> 

POSSIBLE_ALUMINUM 

<|POSSIBLE_AL_PEAKS|>.WAVENUMBER 

441.7 

< 

<|POSSIBLE_AL_PEAKS|>.  WAVENUMBER 

443.7 

> 

«|POSSIBLE_AL_PEAKS|». WAVENUMBER 

470.9 

< 

«|POSSIBLE_AL_PEAKS|».  WAVENUMBER 

472.9 

> 

«<|POSSIBLE_AL_PEAKS|»>. WAVENUMBER 

542.9 

< 

«<|POSSIBLE_AL_PEAKS|>».  WAVENUMBER 

544.9 

Hypothesis:  ALUMINUM_PRESENT 

Actions: 

DO  «|POSSIBLE_AL_PEAKS|».  WAVENUMBER 

AL  CALIBRATION  PEAK.WAVENUMBER 
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created.  In  this  case,  the  hypothesis  ALUMINUM_PRESENT  will  be 
automatically  set  to  false,  and  the  rest  of  the  conditions  will  not  be 
evaluated. 

The  next  six  conditions  use  pattern  matching  to  determine  if 
the  aluminum  confirmation  pattern  is  present.  The  first  two 
conditions  sort  through  the  list  of  objects  in  the  class 
POSSIBLE_AL_PEAKS  to  see  if  a  peak  exists  whose  PEAK_n. WAVE- 
NUMBER  slot  is  greater  than  441.7  and  less  than  443.7.  Since 
aluminum  line  1  has  a  peak  location  of  approximately  442.7 
wavenumbers,  this  narrow  interval  will  be  sufficient  to  locate  it  if 
present.  The  next  two  conditions  look  for  aluminum  line  2,  whose 
wavenumber  is  approximately  471.9.  The  double  set  of  angle 
brackets  tells  the  system  to  pattern  match  on  the  original  list  of 
peaks  in  the  class  POSSIBLE_AL_PEAKS.  If  only  one  set  of  angle 
brackets  was  used,  this  would  mean  "pattern  match  on  the  list 
created  by  the  first  two  conditions,"  which  contains  only  one  peak. 
Since  the  second  two  conditions  look  for  a  completely  different 
peak,  they  need  to  start  over  again  with  the  original  list  of  objects. 
The  last  two  conditions  employ  the  same  methods  to  look  for 
aluminum  line  3,  which  has  a  wavenumber  location  of  543.9 
wavenumbers.  Once  again,  the  triple  set  of  angle  brackets  means  to 
pattern  match  on  the  original  list,  not  the  list  created  by  the  second 
set  of  conditions  with  the  double  set  of  angle  brackets. 

If  all  the  conditions  evaluate  to  true,  then  the  presence  of 
aluminum  is  confirmed  and  the  hypothesis  ALUMINUM_PRESENT  is  set 
to  true.  If  the  hypothesis  is  true,  the  action  will  be  performed.  In 
this  rule,  the  purpose  of  the  action  is  to  assign  the  value  of  the  slot 
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AL_CALlBRATION_PEAK. WAVENUMBER  the  value  contained  in  the 
PEAK_n. WAVENUMBER  slot  that  was  found  with  the  pattern  match 
that  occurred  in  third  and  fourth  conditions  in  the  rule,  which  will 
be  approximately  471.9  wavenumbers.  This  is  indicated  by  the 
double  angle  brackets  around  the  class  |POSSIBLE_AL_PEAKSJ.  The 
'DO'  operator  assigns  the  value  of  the  second  argument  to  the  value 
of  the  first  argument.  So,  after  this  action  is  completed,  the  slot 
AL_CALI BRATION_PEAK.WAV ENUMBER  will  contain  a  value  of 
approximately  471.9  wavenumbers. 

This  calibration  peak  is  used  to  calculate  the  concentration  of 
the  impurity.  Once  the  area  of  the  calibration  peak  is  known,  it  can 
be  multiplied  by  a  calibration  factor  to  determine  the  concentration. 
Each  impurity  has  its  own  calibration  peak,  and  aluminum's  happens 
to  be  the  peak  that  occurs  at  approximately  471.9  wavenumbers.  So, 
after  the  exact  value  of  the  AL_CALIBRATION_PEAK. WAVENUMBER 
slot  is  obtained,  it  is  passed  to  the  FORTRAN  program  'FINDAREA' 
(see  Appendix).  This  program  uses  this  wavenumber  to  locate  the 
peak  in  the  raw  spectral  data  that  is  stored  in  a  common  block.  Once 
it  locates  the  peak,  it  uses  the  trapezoidal  integration  method  to 
determine  its  area.  The  area  is  then  multiplied  by  the  aluminum 
calibration  factor  to  determine  the  aluminum  concentration  in  the 
silicon  sample. 

If  the  aluminum  confirmation  pattern  has  been  identified,  the 
next  step  is  to  identify  all  other  lines  in  the  spectrum  that  are 
caused  by  the  presence  of  aluminum.  Aluminum  may  cause  as  many 
as  20  different  peaks  to  appear,  so  20  different  rules  are  needed  for 
their  confirmation.  Since  the  wavenumbers  where  these  peaks  occur 


is  already  known,  pattern  matching  will  be  used  on  the  class 
|POSSIBLE_ALUMINUM_PEAKS|  to  see  if  they  exist.  This  same 
technique  that  was  used  to  identify  the  aluminum  confirmation 
pattern  will  now  be  used  to  identify  the  individual  lines. 

The  first  rule  will  be  called  ALUMINUM_LINE_1 ,  and  it  will 
attempt  to  confirm  the  presence  of  the  first  aluminum  line,  which 
has  a  peak  location  of  approximately  442.7  wavenumbers.  The  first 
thing  this  rule  must  do  is  check  to  see  if  the  presence  of  aluminum 
has  been  confirmed.  This  is  done  by  examining  the  value  of  the 
hypothesis  ALUMINUM_PRESENT,  which  the  previous  rule  set  to  true 
if  the  aluminum  confirmation  pattern  was  found.  If  ALUMINUM_ 
PRESENT  is  true,  then  the  rule  must  examine  all  peaks  in  the  class 
|POSSIBLE_AL_PEAKS|  to  see  if  one  has  a  peak  location  of  442.7 
wavenumbers.  If  so,  the  slot  AL_LINE_1  .EXIST  will  be  set  to  true 
and  the  slot  AL_LINE_1.WN  will  be  set  to  442.7. 

Name:  Aluminum_Line_1 

Conditions: 

Yes  ALUMINUM_PRESENT 

>  <|POSSIBLE_ALUMINUM_PEAKS|>.  WAVENUMBER  442 

<  <|POSSIBLE_ALUMINUM_PEAKS|>.  WAVENUMBER  443.5 

Hypothesis:  ALUMINUMJJNE 

Actions: 

Let  AL_LINE_1. EXIST  TRUE 

Do  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER 

AL  LINE  1.WN 
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Twenty  rules  like  this  will  exist  for  each  possible  aluminum 
line  that  may  be  present.  For  example,  aluminum__line_X3  has  a  peak 
location  of  approximately  404  wavenumber.  The  rule  to  find  this 
peak  is: 

Name:  AluminumJ_ine_X3 

Conditions: 

Yes  ALUMINUM_PRESENT 

>  <|POSSIBLE_ALUMINUM_PEAKS|>.  WAVENUMBER  403.5 

<  <|POSSIBLE_ALUMINUM_PEAKSl>.WAVENUMBER  404.5 

Hypothesis:  ALUMINUMJJNE 

Actions: 

Let  AL_LINE_X3. EXIST  TRUE 

Do  <|POSSIBLE_ALUMINUM_PEAKS|>.  WAVENUMBER 

AL  LINE  X3.WN 
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3.3 


PROGRAM  FLOW 


The  previous  chapters  demonstrated  how  rules  are  structured 
and  how  they  interact.  Now  it  can  be  shown  how  these  rules  fit  into 
the  flow  of  the  total  program. 

From  now  on,  the  rules  will  be  shown  in  the  way  they  appear  in 
the  rule  network,  which  is  a  graphic  display  of  the  interrelation  of 
all  the  rules.  Looking  at  this  'rule  tree'  allows  the  developer  to 
analyze  the  flow  of  the  program.  The  previous  rules  were  shown  in 
the  context  of  the  rule  editor,  which  is  where  they  are  created. 

For  example,  the  first  rule  the  expert  system  evaluates  is 
named  START_THINGS_OFF: 


Yes  PEAKS_LOADED 
Yes  PEAKS_EXAMINED 
Yes  PEAKS  ANALYZED 


START  THINGS  OFF 


ANALYZED 


This  rule  has  three  conditions  on  the  left-hand  side  that  must 
evaluate  to  TRUE  before  the  hypothesis  ANALYZED  can  be  set  to  TRUE. 
These  three  conditions  represent  three  separate  branches  of  the  rule 
network.  The  purpose  of  this  rule  is  to  force  the  evaluation  of  these 
three  different  branches. 

Th«  first  branch  is  represented  by  the  hypothesis 
PEAKS_LOADED  (Figure  9).  This  branch  of  the  network  loads  the 
spectrum  and  the  peak  table.  After  the  spectrum  and  peak  table  have 
been  loaded,  the  second  branch  of  the  network,  represented  by  the 
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hypothesis  PEAKS_EXAMINED,  is  put  into  action  (Figure  10).  This 
branch  of  the  network  looks  at  the  peak  table  and  classifies  the 
peaks  into  what  they  might  possibly  be,  i.e.,  aluminum,  boron, 
gallium,  indium,  carbon,  oxygen,  or  unknown.  The  next  branch  of  the 
network  is  represented  by  the  hypothesis  PEAKS_ANALYZED  (Figure 
11).  This  branch  takes  the  categories  of  poss’ble  peaks  that  were 
just  created  and  looks  for  the  confirming  patterns  that  tell  which 
impurities  are  present.  The  following  pages  illustrate  how  these 
three  branches  are  connected  to  the  rule  START_THINGS_OFF. 

The  three  branches  of  the  rule  network  interact  in  the 
following  manner.  The  three  conditions  on  the  left  hand  side  point  to 
the  rule  name,  START_THINGS_OFF,  and  the  rule  name  points  to  the 
hypothesis  ANALYZED.  To  get  the  process  going,  the  hypothesis 
ANALYZED  must  be  suggested.  Suggesting  the  hypothesis  forces  the 
expert  system  to  determine  if  its  Boolean  value  is  TRUE  or  FALSE.  To 
do  this,  it  must  evaluate  the  conditions  on  the  left-hand  side. 
ANALYZED  will  be  set  to  true  if  all  three  conditions  evaluate  to  true. 
If  any  of  the  conditions  evaluate  to  false,  ANALYZED  will  be  set  to 
false. 

In  this  rule,  the  left-hand  side  conditions  involve  three 
hypotheses  to  other  rules.  PEAKS_LOADED,  PEAKS_EXAMINED,  and 
PEAKS_ANALYZED  are  all  hypotheses,  and  the  Yes  operator  forces  the 
evaluation  of  their  Boolean  value.  When  the  expert  system  attempts 
to  evaluate  a  hypothesis  that  appears  in  a  condition,  an  event 
called  backward  chaining  takes  place  which  causes  other  rules  in 
that  branch  of  t  he  rule  network  to  be  evaluated.  If  the  value  of 
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Figure  1 1 .  The  PEAKS_ANALYZED  branch  of  the  rule  network. 


the  hypothesis  is  unknown,  then  the  rule  leading  to  that  hypothesis 
must  be  evaluated. 

Since  all  the  rules  in  a  branch  are  connected  and  dependent  on 
each  other,  trying  to  evaluate  one  rule  will  cause  all  rules  leading  up 
to  it  to  also  be  evaluated.  For  example,  the  first  condition  in  the 
rule  START_THINGS_OFF  attempts  to  determine  the  value  of  the 
hypothesis  PEAKS_LOADED,  which  doesn't  yet  have  a  value.  To 
determine  its  value,  the  expert  sy~*em  must  evaluate  the  rule 
LOAD_PEAK_TABLE,  since  this  rule  has  PEAKS_LOADED  as  its 
hypothesis.  If  all  the  conditions  on  the  left-hand  side  of  the  rule 
LOAD_PEAK_TABLE  evaluate  to  true,  then  the  hypothesis 
PEAKS_LOADED  will  be  set  to  true. 

However,  the  first  condition  of  the  rule  LOAD_PEAK_TABLE 
attempts  to  evaluate  the  value  of  SPECTRA_ID_  ASKED  which  is  the 
hypothesis  of  the  rule  GET_SPECTRA_ID.  So  the  rule 
GET_SPECTRA_ID  must  be  evaluated  before  the  hypothesis 
SPECTRA_ID_ASKED  will  have  a  value.  Since  none  of  the  conditions 
in  the  rule  GET_SPECTRA_ID  attempt  to  evaluate  the  hypothesis  of 
another  rule,  this  is  the  end  of  the  'backward  chain.' 

To  summarize,  'suggesting'  the  hypothesis  ANALYZED  forced  the 
evaluation  of  the  rule  START_THINGS_OFF.  The  first  condition  in 
this  rule,  'Yes  PEAKS_LOADED',  forced  the  evaluation  of  the  rule 
LOAD_PEAK_TABLE.  The  first  condition  of  this  rule,  'Yes 
SPECTRA_ID_ASKED',  forced  the  evaluation  of  the  rule 
GET_SPECTRA_ID.  Backward  chaining  forced  the  evaluation  of  all  the 
rules  in  the  PEAKS  LOADED  branch  of  the  rule  network. 
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After  backward  chaining  forces  the  evaluation  of  the  three 
branches  of  the  rule  network,  ANALYZED  will  then  be  set  to  true.  At 
this  point,  the  peak  table  has  been  loaded.  It  has  been  examined  so 
that  the  peaks  are  classified  as  to  what  impurity  they  might  possibly 
represent,  and  the  groups  of  possible  peaks  have  been  analyzed  to  see 
if  the  confirmation  pattern  for  the  various  impurities  exist. 

If  it  has  been  determined  that  the  confirmation  pattern  exists 
for  one  or  more  impurities,  the  next  step  is  to  identify  all  the  peaks 
that  exist  due  to  their  presence.  This  is  accomplished  with  six  more 
rules,  which  are  named  AL_LINES,  BO_LINES,  GA_LINES,  IN_LINES, 
CA_LINES,  and  OXJJNES.  For  example,  if  aluminum  is  present,  then 
the  rule  ALJJNES  will  be  evaluated  to  identify  all  of  the  aluminum 
peaks.  This  rule  will  not  be  evaluated  if  aluminum  is  not  present. 
There  are  approximately  20  different  peaks  that  can  appear  for  each 
impurity,  and  these  rules  will  identify  any  of  these  peaks  that  may 
be  present. 

These  rules  are  evaluated  with  a  process  called  forward 
chaining,  which  is  the  opposite  of  backward  chaining.  In  backward 
chaining,  when  a  hypothesis  appears  as  data  in  a  left-hand  side 
condition  of  some  initial  rule,  the  expert  system  goes  backward 
through  the  rule  network  to  evaluate  all  rules  that  are  needed  to 
determine  the  value  of  this  hypothesis.  In  forward  chaining,  once 
the  value  of  a  hypothesis  or  any  other  data  value  is  determined,  this 
value  is  volunteered  forward  to  any  other  rules  that  uses  it.  Any 
time  the  expert  system  gives  a  rule  a  data  value  it  needs  to  evaluate, 
the  rule  will  automatically  get  evaluated  without  having  to  suggest 
its  hypothesis. 
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Figure  12  illustrates  the  concept  of  forward  chaining.  The 
previous  example  showed  how  three  different  branches  of  the  rule 
network  branched  backwards  from  the  rule  START_THINGS_OFF.  This 
example  shows  how  another  section  of  the  rule  network  branches 
forward  from  the  same  rule.  The  six  rules  that  identify  the  peaks  of 
the  six  impurities  depend  on  the  value  of  the  hypothesis  ANALYZED  of 
the  rule  START_THINGS_  OFF.  This  illustration  shows  how  these  six 
rules  branch  forward  from  the  hypothesis  ANALYZED. 

After  the  rule  START_THINGS_OFF  has  been  successfully 
evaluated,  the  value  for  ANALYZED  will  be  TRUE.  This  value  will  then 
be  volunteered  to  the  six  rules  AL_LINES,  BO_LINES,  GA_LINES, 
INJJNES,  CA_LINES,  and  OXJJNES.  If,  for  example,  the  presence  of 
aluminun  has  been  confirmed,  the  value  of  the  hypothesis 
ALUMINUM_PRESENT  will  be  TRUE.  ALUMINUM_PRESENT  is  the 
hypothesis  of  the  rule  INVESTIGATE_ALUMINUM,  which  looked  for 
presence  of  the  aluminum  confirmation  pattern.  So,  looking  at  the 
rule  AL_LINES,  the  first  condition  'Yes  ANALYZED'  will  be  true,  the 
second  condition  'Yes  ALUMINUM_PRESENT’  will  be  true,  therefore  the 
hypothesis  FIND_AL_LINES  will  evaluate  to  TRUE,  so  the  action 
statement  'Do  ALUMINUMJJNE  ALUMINUMJJNE  will  be  performed. 

ALUMINUM_LINE  is  a  hypothesis,  so  the  action  statement  'Do 
ALUMINUM_L!NE  ALUMINUMJJNE'  is  the  same  as  'suggesting'  this 
hypothesis.  This  simply  forces  the  evaluation  of  any  rule  that 
has  ALUMINUMJJNE  as  its  hypothesis.  In  Figure  13  there  are  twenty 
rules  that  have  ALUMINUM_LINE  as  its  hypothesis. 
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Figure  12.  Forward  chaining  from  the  hypothesis  ANALYZED. 


Each  one  of  these  rules  tries  to  identify  a  particular  peak  and  sets 
the  value  of  the  slot  AL_LINE_/7. EXIST  to  TRUE  if  that  peak  is  found. 
It  also  sets  the  value  of  AL_LINE_n.WN  to  the  wavenumber  of  the 
peak  location.  So,  the  action  statement  'Do  ALUMINUM_LINE 
ALUMINUM_LINE'  triggers  the  evaluation  of  20  rules  which  will 
positively  identify  all  of  the  peaks  that  are  in  the  spectrum  due  to 
the  presence  of  aluminum. 

Following  the  rules  that  identify  the  aluminum  lines  are  the 
rules  that  identify  the  lines  for  all  the  other  impurities  (Figures  14 
through  17).  It  is  not  possible  to  show  all  of  the  rules  for  each  of 
the  impurities.  For  the  rules  that  are  not  shown  in  their  entirety, 
their  peak  locations  can  be  found  in  Table  2. 
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Y*»ALUMtNUM_PRES&rr 
<(POSS _AL,PEAKS)>.WN  »  442.5 
<|POS8jU._PEAKS(>.WN  <  443.5 
•»L*t  AL_LtNE_1. EXIST  TRUE 
-»0o  <)POSS,AL_PEAKS(>.  WN  AL_UNE_1.WN‘ 


Ym  ALUMINUM_PR£SENT 
<(POSS,AL_PEAKSf».WN  >  471.5 
<|POSS_AL_PEAKS(>.WN  <  472.5 
-»Ut  AL_LINE_2.EXIST  TRUE 
-»Do<|POSSJU_PE4KSt>.WN  AL_LINE_2.WN 


Ym  AUJMINUM_PRESWT 
<|PG33_AL_PEAKS|>.WN  >  403.5 
<|POSS _/U._PEAKS|».WN  <  404.5 
->L»t  AL_LINE_X3.  EXIST  TRUE 
>Do<|POSS _AL_PEAKS|>.WN  AL,LJNE_X3.WN 


Ym  AUJMINUMPRESENT 
<|POSS _AL_PEAKS|>.WN  >  413.5 
<|POS3 _AL,PEAKS|».WN  «  41 4.5 
•>L*I  AL_L(NE_X4. EXIST  TRUE 
«.Oo<|POS8J»U_PEAK9|».WN  AL_LINE_X4.WN 

YM  ALUMINUM_PRESENT 
<|POSS,AL_PEAKS|»WN  »  453.5 
<|P03S _AL_PEAKS)>.WN  <  454.5 
-»UPAL_LINE_EX.EXIST  TRUE 
->0o  <|POSS_AL,PEAK8t>.WN  AL_LINE_EX.WN 


ALUMINUMJJNEJ 


ALUMINUM_LINE_2 


ALUMINUM  JJNE_4 
ALUMINUM_LINE_5 


ALUMINUM  JJNE_6 
ALUMINUM  JJNE  7 


ALUMINUM_LINE_8 

ALUMINUM_LINE_9 


ALUMINUM  JJNE_10 
ALUMINUM  _UNE_11 


ALUMINUM_LINE_EI 

ALUMINUM_UNE_2P 


ALUMINUM  LINE  3P 


■  ALUMINUM_LINE_4P 

•  ALUMINUM_UN£_  ESTAR 

■  ALUMINUM JJNE_X1  i 

■  ALUMINUM_LINE_X2  / 


ALUMINUV_LINE_X3 


>ALUMINUM_LINE_X4 


ALUMINUM_LINE_EX 


ALUMINUM  J.INE 


Figure  13.  The  aluminum  rules. 
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Figure  14.  The  boron  rules. 
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Y«»  GALUJM_PRESENT 
<jPOSS_GA_P£AKSf».WN  »  46025 
<|POSS_GA_PEAKS)».WN  <  47055 
”>L«»  GA_UNE_1  .EX ST  TRUE 
-»Oo  «(POSS_QA_PEAKSt».WN  GAJ.NE_I.WN 


V»»  GALLUM_PRESEWT 
<jPOSS_GA_PEAKS(».WN  >  541.0 
<jPOSS_GA_PEAKSl>.WN  <  542.0 
->L*GA_L1NE_3.EXIST  TRUE 
-»Oo  <|POSS_GA_PEAKS)>.WN  GA_LNE_3.WN 


Y««  GALL  UM_PRE  SENT 
<|POSS_GA_PEAKS|>.WN  >  380.0 
<|POSS_GA_PEAKS|>.WN  <  382.0 
-»L«I  GA_UNE_X2.  EXIST  TRUE 
->Oo  <|POSS  .GA_PEAKS)>.WN  GA_UNE_X2.WN. 


Y*»  GALLUMPRESENT 
<|POSS_GA_PEAKS|>.WN  »  422.0 
<|POSS_GA_PEAKS|».WN  «  424.0 
-»L*GA_UNE_X4.EXIST  TRUE 
-»Oo<|POSS_GA_PEAKS|>.WN  GA.LINE_X4.WN 

Y«t  GALLIUM_PRESENT 
<|POSS_GA_PEAKS)>.WN  »  450.0 
<|POSS_GA_PEAKS|>,WN  <  461 .0 
-»L«GAJJNE_EX.EXIST  TRUE 
>»Oo  <|POSS_GA_PEAKSt>.WN  GAJJNE_EX.WN 


GALLIUM  LINE  1 


GALUUM_L!NE_3 


■  GALLIUM_UNE_4  . 

•  GALLIUM_UNE_5  \ 

•  GALLIUM_LIN£_6 

■  GALLIUM_UNE_7 

■  GALLIUM_LINE_8 

■  GALLIUM_LINE_0 

•  GALLIUM  J.INE_10 
.  GALLIUM_LINE_1 1 

■  GALLIUM_LINE_EI 

■  GALUUM_LINE_2P 

•  GALLIUM_UNE_3P 
'  GALLIUM_LINE_4P 
'  GALLIUM_UNE_SP 

■  GALLIUM_LINE_ESTAR 

'  GALUUM_LINE_X1  / 


GALLlUMLIN  E_X2 


GALUUM_UNE_X4 


GALUUM_LINE  EX 


GALLIUM  LINE 


Figure  15.  The  gallium  rules. 
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Yes  INDIUM_PRESENT 
<|POSSJN_PEAKS|>.WN  >  1 145.0 
<]POSS_IN„PEAKS|>.WN  <  1 146.0 
->Let  IN_LINE_1. EXIST  TRUE 
>Do  <|POSS_IN_PEAKS|>.WN  IN_LINE_1  .WN 


Yes  INDIUM_PRESENT 
<|POSSJN_PEAKS|>.WN  >  1 175.3 
<|POSS_IN_PEAKSj>.WN  <  1 176.3 
->Le(  INJJNE_2.EXIST  TRUE 
->Do  <lPOSS__IN_PEAKS|>.WN  IN_LINE_2.WN 


Yes  lNDIUM_PRESENT 
<|POSS_IN_PEAKS|>.WN  >  862.0 
<|POSS_IN_  PEAKS!  >.WN  <  864.0 
->let  !N_LINE_X3. EXIST  TRUE 
-»Do  <|POSS_IN_PEAKS|>.WN  IN_UNE_X3.WN 


Yes  1NDIUM_PRESENT 
<|  POSS  JN_P£AKS|>.WN  >  872.0 
<|POSS_lN_PEAKS|>.WN  <  874.0 
->Let  IN_UNE_X4. EXIST  TRUE 
->Do  <|POSS_IN_PEAKS|>.WN  IN_LINE_X4.WN 


Yes  IN0IUM_PRESENTT 
<|POSS_IN  PEAKS|>.WN  >  908.0 
<|POSS  _IN  PEAKS|>.WN<  910.0 
->L«HN_UNE_ EX. EXIST  TRUE 
Do  <|POSS_IN_PEAKS|>.WN  1NJ.INE_EX.WN 


INDIUM  JJNE_1 


INDlUM_UNE_2 


INDIUM  JJNE_3 
INDIUM_LIN6_4 


iND>UM_LINE_4A 

INDIUM_L1NE_5 


!NDIUM_LINE_6 

INDIUM_UNE_7 


INDIUM_LINE_8 
INDIUM  LtNE_9 


1NDIUM_LINE_10 
INDJUM_LINE_1 1 


IND!UM_LINE_EI 
INDIUM  LINE  2P 


INOIUM_LINE_3P 


INDIUM_LINE_4P 
'  IND!UM_L!NE_ ESTAR 
•  INDlUM_UNE_Xt 


INDIUM  UNE_X2 


INDIUM.  UNE_X3 


IN01UM_LINE_X4 


IN0lUM_LINE_EX 


INDIUM  LINE 


Figure  16.  The  indium  rules. 


CARSON  LINE 


Figure  17.  The  carbon  and  oxygen  rules. 
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4.0  SAMPLE  RUN 


Table  3  is  an  example  of  the  report  that  is  produced  by  the 
expert  system.  For  this  sample  run,  it  was  given  a  silicon  spectrum 
that  contained  the  impurities  boron,  gallium,  and  oxygen.  As  seen  in 
the  report,  the  expert  system  was  able  to  identify  these  impurities, 
determine  their  concentrations,  and  list  their  lines  that  were 
present. 
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Table  3.  The  report  generated  by  the  expert  system. 


RESULTS  OF  IMPURITY  CLASSIFICATION 

11/13/1990 

SPECTRUM  NAME:  siga819ca 
TYPE:  abs 


IMPURITY: 

BORON 

CONCENTRATION: 

2.72301772882080081615 

PEAKS  DETECTED 

LINE 

WAVENUMBER 

1 

244.96 

2 

278.23 

2p 

669.0 

3 

309.5 

4 

320.08 

8 

346.33 

IMPURITY: 

GALLIUM 

CONCENTRATION: 

3.8677543243408203e1 6 

PEAKS  DETECTED 

LINE 

WAVENUMBER 

1 

469.85 

2p 

897.04 

3 

541.96 

3p 

921.13 

4 

550.03 

6 

569.19 

8 

574.21 

10 

582.28 

XI 

346.33 

IMPURITY: 

OXYGEN 

CONCENTRATION: 

2.581 75025939941 4050e1 6 

PEAKS  DETECTED 

LINE 

WAVENUMBER 

1 


1136.92 
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5.0  TESTING  THE  EXPERT  SYSTEM 


In  deciding  what  should  be  considered  in  testing  the  expert 
system,  four  main  criteria  were  settled  upon: 

1.  Easy  to  use:  The  system  should  be  easy  to  use  with  a  minimum 
of  experience  needed  in  computer  use. 

2.  Requires  little  analysis  experience:  The  expert  system  should 
produce  satisfactory  results  when  being  used  by  a  lower  level 
researcher  such  as  a  technician,  co-op  student  or  entry-level 
researcher. 

3.  Produces  accurate  results:  The  results  generated  by  the  expert 
system  should  agree  with  those  of  an  experienced  spectroscopist. 
These  results  are  often  quite  subjective  and  vary  considerably  when 
the  analysis  is  performed  by  different  researchers.  The  expert 
system  should  find  all  the  impurities  present  in  a  given  sample 
spectra  and  the  calculated  impurity  concentrations  should  be  within 
the  experimental  tolerances  that  are  considered  acceptable  in  this 
type  of  spectral  analysis. 

4.  Produces  consistent  results:  The  expert  system  should  give 
consistent  results  if  given  identical  spectra.  This  may  seem  like  an 
obvious  requirement,  but  the  user  is  allowed  some  interaction  with 
the  system  and  can  add  peaks  to  the  peak  table  if  he  feels  they  are 
significant  when  viewing  the  spectrum. 
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5.1  DISCUSSION 


The  expert  system  was  first  tested  by  an  experienced 
spectroscopist  with  little  computer  experience.  No  significant 
difficulties  were  experienced  by  this  user.  After  one  tutorial 
session,  the  researcher  was  confident  enough  to  run  the  system 
without  assistance.  A  co-op  student  was  the  next  test  subject  and 
required  only  one  short  tutorial  session  to  become  confident  enough 
to  use  the  system  without  assistance.  The  only  interaction  required 
by  the  expert  system  is  for  the  user  to  s'  pply  a  spectrum  name  and 
indicate  whether  it  is  an  absorption  coefficient  or  photoconductivity 
spectrum.  Then  the  user  is  asked  to  examine  the  spectrum  on  the 
screen.  This  spectrum  has  all  peaks  marked  that  the  expert  system 
finds.  The  user  may  mark  additional  peaks  if  he  feels  they  are 
significant.  If  these  peaks  are  simply  noise  or  some  other 
irregularity  in  the  spectrum,  then  the  expert  system  ignores  them. 
This  is  all  the  user  is  required  to  input  to  the  system. 

In  order  to  test  the  accuracy  of  the  expert  system,  its  results 
were  compared  with  the  results  of  a  study  conducted  several  years 
earlier.  This  study  examined  the  relationship  between  the  area  under 
a  specific  gallium  peak  in  the  spectrum  and  the  gallium 
concentration  for  that  sample  as  determined  by  other  means  such  as 
Hall  Analysis  or  Photoluminesence.  The  researcher  conducting  this 
study  ran  a  series  of  gallium  doped  samples  and  carefully  measured 
what  he  felt  the  gallium  dopant-concentration  peak  areas  were.  The 
results  of  this  study  were  used  to  calculate  a  concentration 
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coefficient  relating  a  specific  gallium  peak’s  area  to  the  gallium 
concentration  present  in  the  sample. 

Fourteen  spectra  were  analyzed  by  the  expert  system  in  this 
test.  Some  runs  were  made  by  us,  the  designers,  others  by  a  co-op 
student  working  in  the  lab.  The  results  are  displayed  in  table  4. 


Table  4.  Results  of  gallium  concentration  study  test. 


SAMPLE  ID 

GALLIUM  CONCENTRA 

GA0102AC 

-4.33 

GA044BAC 

-5.66 

GA1004AC 

-3.23 

GA1069BAC 

-1 .46 

GA1305AAC 

-7.0 

GA203AC 

+0.6 

GA370AC 

+2.1 

GA402AC 

-6.6 

GA757DAC 

+  0.1 

GA757EAC 

-2.5 

GA757FAC 

+0.42 

GA987AC 

+1.25 

GA987BAC 

-8.96 

GA1072AC 

-5.26 

The  expert  system's  results  for  impurity  or  dopant 
concentration  averaged  2.9%  less  than  the  values  calculated  by  the 
spectroscopist  in  this  study.  This  figure  is  well  within  a  reasonable 
value  if  comparing  the  results  of  two  experienced  researchers. 
Typically  an  error  of  20%  is  not  unreasonable  due  to  differences  that 
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occur  when  determining  the  area  of  the  calibration  peak.  A 
calibration  peak  is  a  peak  occurring  in  a  spectrum  that  is  always 
generated  by  the  presence  of  a  particular  dopant  or  impurity.  The 
area  of  this  peak  is  multiplied  by  a  scale  factor  to  calculate  that 
particular  dopant's  concentration.  If  the  individual  chooses 
different  peak  boundaries,  then  different  impurity  or  dopant 
concentrations  will  result.  In  a  noisy  spectrum,  peak  boundaries  can 
be  difficult  to  discern.  This  situation  also  occurs  when  the 
concentrations  are  very  low  and  the  peak  is  relatively  weak.  This 
can  be  seen  in  Figure  18;  a  typical  weak  peak  occurring  in  a  noisy 
area  of  the  spectrum.  Defining  the  limits  of  a  peak  in  a  noisy  area 
can  be  difficult  at  best.  This  is  where  the  differences  generally 
occur  when  two  researchers  analyze  the  same  spectrum. 

The  dopant  or  impurity  peaks  were  classified  correctly  in  all 
cases.  The  expert  system  always  detected  the  presence  of  all 
impurities  known  to  be  present  in  the  samples.  In  summary,  there 
was  good  agreement  between  the  system's  analysis  and  those  of 
other  experienced  researchers. 
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BORON 


260  270  280  290 

WAVENUMBER 

Figure  18.  Possible  peak  boundaries  in  a  noisy  spectrum. 


55 


6.0  FUTURE  PLANS  FOR  THE  EXPERT  SYSTEM 


The  expert-system  is  completed  for  use  with  high  purity 
silicon  materials.  As  stated  in  previous  sections,  the  scope  of  the 
project  was  limited  to  silicon  based  electronic  materials.  The  body 
of  knowledge  for  these  materials  was  well  established.  To  have 
included  other  types  of  optical  materials  in  the  knowledge  base 
would  have  required  more  man-hours  of  development  and  testing 
than  could  be  justified.  The  main  goal  was  to  use  this  as  a  learning 
experience  for  future  expert-system  work  within  our  laboratory. 
The  completed  system  is  dependable  and  could  be  used  on  a  day-to- 
day  basis  for  research  as  well  as  quality-control  applications. 
Unfortunately,  silicon  based  electronic  materials  are  no  longer  the 
current  materials  of  interest,  so  it  is  unlikely  the  system  will  be 
heavily  used. 

The  next  phase  of  this  project  is  to  design  a  knowledge  base 
for  analyzing  superlattice  materials  grown  by  Molecular  Beam 
Epitaxy  (MBE).  There  currently  is  not  a  body  of  knowledge  to  draw 
upon  for  constructing  this  new  system.  When  a  sufficient  number  of 
samples  of  this  type  have  been  analyzed,  a  set  of  rules  should 
emerge  which  will  be  the  basis  for  the  next  phase.  There  has  been 
discussion  of  designing  a  closed-loop  system  where  the  expert 
system  can  control  the  spectrometer  as  well  as  the  sample  chamber 
temperature,  but  this  plan  is  stiil  tentative  and  requires  much  more 
study  as  well  as  the  full  support  of  the  researchers  in  this  area. 


56 


APPENDIX  A.  FORTRAN  SUBROUTINES 


The  expert-system  calls  several  FORTRAN  subroutines  to 
perform  routine  analysis  tasks.  These  routines  are  described  below  . 

FILLSPECTRA  The  spectral  array  of  data  points  is  stored  in  a  named 
Fortran  common  block.  FILLSPECTRA  is  passed  the  name  of  the 
spectra  to  be  analyzed.  It  simply  locates  the  the  spectra  data  file, 
opens  it  and  loads  the  common  block  with  the  data  points  from  this 
file.  These  data  points  are  subsequently  used  for  the  interractive 
graphics  in  PEAKPICK  and  in  calculations  and  graphics  in  FINDAREA. 

PEAKPICK  This  is  an  interractive  graphics  routine  which  plots  the 
subject  spectrum  on  the  workstation  screen.  All  peaks  that  have 
been  located  by  the  FTIR  data  station  are  marked  and  the  user  is  then 
asked  if  he  wishes  to  mark  additional  peaks  of  interest.  By  using  the 
attached  mouse  and  cursor  all  peaks  of  interest  that  were  not 
marked  by  the  program  can  be  marked  by  the  user  interractively. 
When  complete  this  routine  transfers  to  a  subroutine  called  REFORM 
which  creates  an  accompanying  table  of  peaks  from  the  spectrum 
under  scrutiny.  This  table  is  in  a  format  readable  by  the  expert- 
system  program.  When  the  knowledge  base  is  loaded  this  peak  table 
is  used  as  a  basis  for  generating  preliminary  conclusions  about  the 
sample. 
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F1NDAREA  The  impurity  concentration  is  calculated  by  multiplying 
the  area  of  a  certain  peak  associated  with  that  impurity  by  a 
calibration  factor.  The  resulting  impurity  concentration  is  in  the 
units  of  atoms/cm**2.  Since  these  peaks  are  superimposed  on  a 
background  signal  that  is  often  noisy,  or  the  peak  itself  is  very 
weak,  defining  the  peak  limits  can  be  a  subjective  process.  Different 
researchers  can  analyze  the  same  spectra  and  obtain  peak  areas 
differing  by  20%.  This  routine,  after  having  been  passed  the 
wavenumber  of  an  impurity  peak  of  interest,  displays  suggested  peak 
boundaries  and  asks  the  user  if  he  agrees:  if  yes  then  calculations 
continue  based  on  these  peak  limits.  If  no,  the  user  indicates  with 
the  mouse  and  cursor  what  the  peak  boundaries  should  be  and  these 
are  used  for  subsequent  calculations.  The  peak  area  is  then 
calculated  by  a  trapezoidal  integration  method.  The  user  interaction 
is  handled  by  a  subroutine  called  PEAKPLOT,  called  from  FINDAREA. 
The  calculated  value  of  the  area  of  the  peak  of  interest  is  then 
passed  back  to  the  expert-system  to  be  recorded  and  used  in 
impurity  concentration  calculations.  All  FORTRAN  callable  screen 
graphics  are  produced  by  a  graphics  software  package  called  CA- 
DISSPLA,  distributed  by  Computer  Associates  Inc. 
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FILLSPECTRA 


This  routine  reads  in  the  spectrum  points  and  loads  them  into 
a  common  block  called  BLOCK1. 


INTEGER  FUNCTION  FILLSPECTRA(THESTR,NATOMS,THEATOMS) 
INCLUDE  'NXP$LIBRARY:NXPVMSDEF.FOR* 

CHARACTER*!*)  THESTR 
INTEGER*4  THEATOMS(*) 

INTEGER*4  RETVAL 

LOGICAL  FLAG  I  if  returned  TRUE  then  no  errors 
INTEGER’4  NATOMS 
CHARACTER  ANS 
CHARACTER*80  LINE 

CHARACTER*10  INFILE  IThe  datafile  with  no  extension 
CHARACTER*^  DATAFILE  I  the  spectra  file 

COMMON  /BLOCK  1  /X(20000) ,  Y(20000),NUMBPTS,INFILE,FLAG 
I  the  spectral  array  of  X-Y  pairs 


RETVAL=NXP$GETATOMINFO(THEATOMS(1),NXP$_AINFO_VALUE,0,0, 
1  .INFILE, NUMB1  )  I  the  spectra  file  name 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

open  the  file  to  be  examined 

This  section  reads  in  the  spectral  array 


c 

DATAFILE  =  INFILE//\MAS' 

OPEN  (FILE=DATAFILE,STATUS=OLD'1UNIT=1,CARRIAGECONTROL='UST', 
1  ERR=5000) 

C  the  following  section  reads  the  input  file's  header  section 
C  if  there  is  one 
C 

READ(1 ,30)  LINE 
3  0  FORMAT(A80) 

DO  WHILE  (LINE(1 :2).EQ.'##') 

READ(1 ,30)  LINE 
END  DO 

get  the  actual  data  X-Y  pairs 
1=1 


C 


80 


READ{1  END-95)  X(I),Y(I) 

1-1+1 

NUMBPTS  -  I 

GO  TO  80  I  read  until  an  EOF  is  found 
9  5  CLOSE(  UNIT-1  )  I  normal  exit  point 
FILLSPECTRA  -  1  I  return  success 
RETURN 

5000  print  5001 

5001  formate  BAD  FILE  NAME ') 
return 
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PEAKPICK 


INTEGER  FUNCTION  PEAKPICKQ 


This  routine  allows  the  user  to  view  the  spectrum  with  the 
peaks  marked.  The  user  can  then  interactively  mark  add¬ 
itional  peaks  that  he  feels  should  be  included  in  the 
analysis  of  the  sample. 


***************  Variable  Declaration  * . ************* 

X,Y  The  x  &  y  values  of  the  whole  spectra 

Numbpts  The  number  of  points  in  the  spectra 

Num  1/4  of  the  number  of  peaks  in  the  whole  spectra 

Npkl  -4  The  number  of  points  in  each  corresponding  plot 

XI  ,Y1  The  temp  var  with  the  points  to  be  plotted 

for  each  of  the  four  plots 

Peakl  -4  The  peaks  picked  by  the  user  for  each  plot 

Oldpkx.OldpkY  The  peaks  found  in  the  .pk  file 
Name  The  raw  name  of  the  spectra  (no  exts) 

Line  Dummy  variable  used  to  read  in  a  header  line 

INCLUDE  'NXP$L!BRAR  Y  :NXP  VMSD  EF.FOR’ 

Common  /Block1/X(20000),  y(20000),  Numbpts,  Name, flag 

Integer  Num, Npkl  ,Npk2,Npk3,Npk4,numold 

real  XI  (5000),  Y1(5000) 

real  Peakl  (25,2),Peak2(25.2),Peak3(25,2),Peak4(25,2) 

real  OldPkx(IOO),  OldPky(IOO) 

character*  10  Name 

character*15  infile 

character*80  Line 

logical  flag,  exist 


Begin  Code 


IF  (EXIST(name//'.pk'))  THEN 

Open(unit=1 ,  statu  s='old',file=name//'.pk’,err=5000) 
I  Loop  to  get  rid  of  header  in  .pk  file 
do  20  1-1,6 
read(1,10)  line 
format(a80) 
continue 


I  Loop  to  read  in  the  old  peaks  from  the  .pk  file 
1  =  1 

read(1  ,*,end=45)  OldPkx(l),  OldPky(l) 

1=  I  +  1 
go  to  30 
continue 
numold  =  I 


else 


open(unit=1  ,status='new,,{ile=name//‘.pk’) 
do  50  1=1,6 

write(1,*)  'Header  Line' 

5  0  continue 

Numold  =  0 

endif 

C**  Here  we  send  the  plot  to  the  screen  in  four  different  parts,  having 
C**  the  user  return  the  peaks  he  sees  on  each  of  the  four  plots 


Num  =  Numbpts/4 

C******  First  Plot 
K=0 

Do  60  1=1, Num 
K=K+1 
x1(k)=x(l) 

Y1(k)=Y(l) 

6  0  continue 

call  Plot(X1  ,Y1  ,K,Name,Peak1  ,Npk1  .oldpkx.oldpky, numold) 

C******  Second  Plot 
K=0 

Do  70  l=(Num-100),(Num*2) 

K=K+1 

x1(k)=x(l) 

Y 1  (k)=Y  (I) 

7  0  continue 

call  Plot(X1  ,Y1  ,K, Name, Peak2,Npk2, oldpkx.oldpky, numold) 
Third  Plot 

K=0 

Do  80  l=(Num*2-100),(Num*3) 

K=K+1 

x1(k)=x(;) 

V1(k)=Y(l) 

8  0  continue 

call  Plot(X1  ,Y1  ,K, Name, Peak3,Npk3, oldpkx.oldpky, numold) 

C******  Fourth  Plot 

K= 0 

Do  90  l=(Num*3-100),Numbpts 

K=K+1 
xl  (k)=x(J) 

Y 1  (k)= Y  (I) 

9  0  continue 

call  PlotfXI  ,Y1  ,K)Name,Peak4,Npk4Ioldpkx,oldpky, numold) 

call  Donepl  Ifinished  plotting 
CALL  SETDEV(IERR.O) 
close(l)  Iclose  the  .pk  file 
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C****  Now  we  take  the  four  peak  arrays  and  append  them  to  the  Name.pk 
C*“*  file  that  already  exist 

infile  =  name//'.pk' 

IF  (EXIST (infile))  then 

Open(unit=1 0,status=,old,1  access='append',file=infile) 

ELSE 

Open(unit=1 0,status='new',  file=infile) 

END  IF 

do  500  1=1  ,Npk1 

Write(10,*)  peakl(l.l),  peak1(l,2) 

5  00  Continue 

do  510  1=1  ,Npk2 

Write(10,*)  peak2(l,1),  peak2(l,2) 

5 1  0  Continue 

do  520  1=1  ,Npk3 

Write(10,*)  peak3(l,1),  peak3(l,2) 

520  Continue 

do  530  1=1  ,Npk4 

Write(10,*)  peak4(l,1),  peak4(l,2) 

530  Continue 

close(IO)  (close  the  .pk  file 
PEAKPICK  =  1  !  return  success  code 

CALL  REFORM  I  make  a  new  .nxp  peak  table 

return 
c 
c 

5000  wr'*e(V)  name//'.pk',  'error  opening  that  file!!' 

FLAL>  =  .FALSE, 
call  exit 
end 


C 


PLOT  . . . ***** . * 

subroutine  Plot(X,Y,N1Name,Peak,NumPts,oldpkx,oldpky,numold) 

real  X(5000),  Y(5000),  Peak(100),olpkx(100),oldpky(l00) 

real  Xmin,Xstep,Xmax,Ymin,Ystep,Ymax 

Integer  N,  NumPts,  NumOld 

Character'10  Name 


N=N- 1 
call  PGPX 

call  newclr('whit’) 
call  hwspec(0.0,  'BACK') 
call  page(13,11.) 
call  area2d(12,10) 
call  NEWCLR('viol’) 
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call  Headin(name,  100,1. 5,1) 

call  Rndlin(X,N,1 1 .5,Xmin,Xstep,Xm  > <} 

call  Rndlln(Y,N,9.5IYmin,Ystep,Ymax) 

call  Xname('Wave  Number',100) 

call  Yname(’RangeMOO) 

call  Graf(Xmin,Xstep,Xmax,Ymin,Ystep,Ymax) 

call  Newclr('red') 

call  Curve(X,Y,N,0) 

call  N8wclr('Gree’) 

call  Marker(13) 

call  Setdev(0,0) 

call  Curve(oldpkx,oldpky,Numold,-l ) 

Call  messag('Select  peaks  with  left  mouse  button. ’,37,1 .0,9.75) 

Call  Messag(’Click  middle  button  to  continue. ’,31 ,7.5,9.75) 

call  Read_Screen(Xmin,Xmax,Ymin,Ymax,Peak,Numpts) 

call  Endpl{0) 

return 

end 


C 


Subroutine  Exist 


Function  Exist(FileName) 

Tests  for  the  existance  of  a  filename. 

Character*15  FileName 
Integer  IERR 
Logical  Exist 

Test  for  file  at  CurDate 

OPEN  (File  =  FileName, status  =  'OLD’, Unit  =  l.iostat  =  ierr) 
if  (IERR.EQ.29)  then 
Exist  =  .FALSE, 
else 

Exist  =  .TRUE, 
end  if 


Close(Unit  =  1) 

Return 

End 

q. ........... .  Read  screen 


subroutine  read_screen  (xmin,xmax,ymin,ymax,Peak,NumPts) 

real  xpos  lx  position  read  from  the  screen 

real  ypos  ly  position  read  from  the  screen 

integer  i  Inumber  of  points  read  form  the  screen 

real  save_xpos(100)  Ithe  numbers  inputed  from  the  screen 
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real  save_ypos(100)  Ithe  numbers  inputed  from  the  screen 
logical  first  lis  this  the  first  pos  read  from  the  screen? 

integer  rkey  Istatus  of  the  return  key  (0=o.k.,  1=stop) 

real  Peak{25,2) 

Integer  NumPts 

C  Use  the  mouse  to  get  the  positions  from  the  screen 
C  call  angle(9.0) 
first  =  .true, 
rkey  =  0 
1  =  0 

do  while  (rkey  .ne.  1) 

1 0  continue 

call  get_position(xmin,xmax,ymin,ymax, first, 

&  xpos, ypos, rkey) 

if  (rkey  .eq.  0)  then 
call  MARKER(13) 
call  CURVE(Xpos,Ypos,1  ,-1 ) 

1  =  1  +  1 

Peak(l.l)  =  Xpos 
Peak(l,2)  =  Ypos 
save_xpos(l)  =  xpos 
save_ypos(l)  =  ypos 
go  to  10 

else  if  (rkey  .eq.  1)  then 
NumPts  =  I 
return 
else 

type  300 

300  formate  No  valid  data  available') 

endif 
enddo 

call  ANGLE(O.O) 

return 

end 


C*** . ***********  g  exposition 

subroutine  GET_POSITION(xmin,xmax,ymin,ymax, 

&  first.xpos,  ypos,  rkey) 

Integer  rkey  Istatus  return  key  (0=ok,1=stop) 

Real  xmin,xmax,ymin,ymax 

Real  pos,  Ypos  Ireturned  values  from  the  plot 

Logical  First 

Integer  IKEY_R 

Real  Initx,  Inity  {initial  cross-hair  position 


if  (first)  then 
call  CRDSYS('USER') 
initx=(xmax-xmin)/2.0 


inity-(ymax-ymin)/2.0 
firsts  -  .false, 
endif 

call  INILOC(1,xpos,ypos,0,0) 

call  ECHL0C(1 ,2,xmin,xmax,ymin,ymax) 

call  REQL0C(1  ,rkey,xpos,ypos,ikey_r) 

return 

end 
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FINDAREA 


This  routine  is  passed  the  left  and  right  x-coordinates  of  a  peak  in  a  spectral 
data  file.  This  subroutine  calculates  the  area  of  that  peak  and  returns  it  to 
the  expert  system  for  later  use. 


SUBROUTINE  FINDAREA(THESTR,NATOMS,THEATOMS) 

INCLUDE  ,NXP$UBRARY:NXPVMSDEF.FOR’ 

CHARACTER***)  THESTR 
INTEGER*4  THEATOMS(*),STAT 

INTEGER*4  RETVAL.RETVAL1  ,RETVAL2,RETVAL3,IDFLAG 
LOGICAL  FLAG  I  if  returned  TRUE  then  no  errors 
INTEGER‘4  NATOMS 

COMMON  /BLOCK1/X(20000),  Y(20000), LOADED, NUMBPTS 
I  the  spectral  array  of  X-Y  pairs 

CHARACTER*80  LINE 

CHARACTER'1 5  INFILE  !  the  spectra  file 
CHARACTERS  INCODE 
CHARACTER*2  CODE 
REAL  LASTX,MR,ML,L_EDGE,PEAK 
INTEGER  DELTA, RINDEX 


RETVAL1  =NXP$GETATOMINFO(THEATOMS(1  ),NXP$_AINFO_VALUE,0,0, 
$  , INFILE, NUMB1  )  I  the  spectra  file  name 

RETVAL2=NXP$GETATOMINFO(THEATOMS(2),NXP$_AINFO_VALUE,0,Q, 
$  NXP$_DESC_FLOAT,PEAK,NUMB2  )  !  the  peak  location 

RETVAL3=NXP$GETATOMINFO(THEATOMS(3),NXP$_AINFO_VALUE,0,0, 
$  .INCODE, NUMB3  )  !  the  type  of  dopant  code 

IF  (RETVAL1  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 
IF  (RETVAL2  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 
IF  (RETVAL3  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

CODE=  INCODE{1:NUMB3)  I  cut  off  excess  stuff 


IPEAKX=IFIND(PEAK)  !  get  set  up  for  loop 

IRIGHT=  0 

LEFT=  10 

MR  =  -20 

ML=  20 

C  This  loop  works  its  way  down  both  sides  of  the  peak,  computing  the 
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C  average  slope  on  either  side.  When  one  side  bottoms  out  (  slope  is 
C  zero  )  it  stops. 

DO  WHILE(  (ML.GT.O)  .AND.  (MR.LT.O)) 

MR.  SLOPE(IPEAKX+IRIGHT)  I  get  the  right  side  slope 
MU  SLOPE(IPEAKX-LEFT)  I  get  the  left  side  slope 
LEFT=  LEFT+1  I  farther  down  left 

IR1GHT.  1RIGHT+1  I  farther  down  right 
END  DO 

IF  (ML.LE.O)  DELTA.  LEFT  +  6  I  why  6  and  4,  you'll  have  to 
IF  (MR.GE.O)  DELTA.  IRIGHT  +  4  I  draw  a  picture  and  convince 

I  yourself  that  these  work 

LINDEX.  IPEAKX  -  DELTA  I  the  left  edge  index  of  the  peak 
RINDEX.  IPEAKX  +  DELTA  I  the  right  edge  ..idex 
L_EDGE  .  X(LINDEX) 

R_EDGE  .  X(RINDEX) 

C 

CALL  PEAKPLOT(L_EDGE,R_EDGE,CODE)  !  let  the  user  look  at  the  peak 

I  limits  &  change  them  if  needed 

CALL  PKAREA(L_EDGE,R_EDGE,AREA)  I  get  the  area 

IF  {  AREA  .LE.  0)  GO  TO  5000  I  if  its  a  bad  area,  say  so 


IF  (C ODE.EQ.'AL’)  THEN  I  if  it's  an  aluminum  peak 

RETV  AL=NXP$GETAT0MID('ALUMINUM_CAUBRATI0N_PEAK  .AREA', 

1  IDFLAG.NXP$_ATYPE_SLOT)  I  get  the  area  address 
IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

RETVAL=NXP$GETATOMID('ALUMINUM_CAUBRATION_PEAK.START', 
1  IDFLAG1  ,NXP$_ATYPE_SLOT)  I  get  the  l  edge  address 
IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

RETVAUNXP$GETATOMID('ALUMINUM_CALIBRATION_PEAK.END', 

1  IDFLAG2,NXP$_ATr  PE_SLOT)  I  get  the  R-edge  address 

IF  (RETVAL  .NE.  NXP$  _ERR_NOERR)  GO  TO  5000 

c  volunteer  the  area  value 

RETV  AUNXP$VCLUNTEER(IDFLAG,NXP$_DESC_FLOAT, AREA, 

1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 
C 

c  volunteer  the  l-edge  value 

C  RETVAUNXP$VOLUNTEER(!DFLAG1  ,NXP$_DESC_FLOAT,L_EDGE, 

C  1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

C 

C  volunteer  the  r-edge  value 

C  RETVAUNXP$VOLUNTEER(IDFLAG2,NXP$_DESC_FLOAT,R_EDGE, 
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c 

c 


1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 
IF  {RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 
GO  TO  95 
ENDIF 


IF  (CODE.EQ.'CA')  THEN  !  if  it’s  an  carbon  peak 
RETVAL=NXP$GETATOMlD(,CARBON_CAUBRATION_PEAK.AREA\ 

1  IDFLAG,NXP$_ATYPE_SLOT)  !  get  the  area  address 
IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

RETVAL=NXP$GETATOMID(,CARBON_CAUBRATION_PEAK.START', 
1  IDFLAGl , NX P$_AT Y P E_S LOT)  !  get  the  l-edge  address 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

RETVAUNXP$GETATOMID('CARBONCALlBRATION_PEAK.END,1 
1  IDFLAG2,NXP$_ATYPE_SLOT)  I  get  the  R-edge  address 
IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

volunteer  the  area  value 

RETVAL=NXP$VOLUNTEER(IDFLAG,NXP$_DESC_FLOAT,AREA, 

1  JlOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

volunteer  the  l-edge  value 

RETVAL=NXP$V0LUNTEER(IDFLAG1,NXP$_DESC__FL0AT,L_EDGE, 
1  JlOR(NXP$_VSTRAT_NOFWRD.NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

volunteer  the  r-edge  value 

RETVAL=NXP$VOLUNTEER(IDFLAG2,NXP$_DESC_FLOAT,R_EDGE, 
1  JlOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 
GO  TO  95 
ENDIF 


IF  (CODE.EQ.’OX’)  THEN  !  if  it’s  an  oxygen  peak 
RETV  AL=NXP$GETATOMID('OXYGEN_CALIBRATION_PEAK  .AREA', 

1  IDFLAG,NXP$_ATYPE_SLOT)  I  get  the  area  address 
IF  (RETVAL  .NE.  NXP$_ERR_NOERR;  GO  TO  5000 

C  RETVAUNXP$GETATOMID('OXYGEN_CAUBRATION_PEAK.START', 

C  1  IDFLAGl  ,NXP$_ATYPE_S LOT)  I  get  the  l-edge  address 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

C  RETVAL=NXP$GETATOMID('OXYGEN_CALIBRATION_PEAK.END’, 

C  1  IDFLAG2,NXP$_ATYPE_SLOT)  I  get  the  R-edge  address 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

c  volunteer  the  area  value 

RETVAL=NXP$VOHJNTEER(IDFLAG,NXP$_DESC_FLOAT,AREA, 
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1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

c  volunteer  the  l-edge  value 

C  RETVAL-NXP$VOLUNTEER(IDFLAG1  ,NXP$_DESC_FLOAT,L_EDGE, 

C  1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

C  volunteer  the  r-edge  value 

C  RETVAL=NXP$VOLUNTEER(lDFLAG2,NXP$_DESC_FLOAT,R_EDGE, 

C  1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

GO  TO  95 
ENDIF 


IF  (CODE.EQ.'BO')  THEN  I  if  it's  a  boron  peak 
get  the  area  address 

RETVAL«NXP$GETATOMID('BORON_CALIBRATION_PEAK.AREA,1 
1  IDFLAG,NXP$_ATYPE_SLOT) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

RETVAL=NXP$GETATOMID('BORON_CALIBRATION_PEAK.START', 

1  IDFLAG1  ,NXP$_ATYPE_SLOT)  I  get  the  l  edge  address 
IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

RETVAUNXP$GETATOMID('BORON_CAUBRATION_PEAK.END,> 

1  IDFLAG2,NXP$_ATYPE_SLOT)  !  get  the  R-edge  address 
IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

volunteer  the  area  value 

RETVAUNXP$VOLUNTEER(IDFLAG,NXP$_DESC_FLOAT,AREA, 

1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

volunteer  the  l  edge  value 

RETVAUNXP$VOLUNTEER(IDFLAGl,NXP$_DESC_FLOAT.L_EDGE, 
1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 


volunteer  the  r-edge  value 

RETVAUNXP$VOLUNTEER(IDFLAG2,NXP$_DESC_FLOAT,R_EDGE, 
1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 
GO  TO  95 
ENDIF 


IF  (CODE.EQ.'IN')  THEN  !  if  it's  an  indium  peak 
c  get  the  area  address 
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RETVAL=NXP$GETATOMID(,lNDIUM_CALIBRATION_PEAK.AREA', 

1  IDFLAG,NXP$_ATYPE_SLOT) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

C  RETVAL=NXP$GETATOMID('INDIUM_CALIBRATION_PEAK.START’, 

C  1  IDFLAG1  ,NXP$_ATY P E_S LOT)  I  get  the  l-edge  address 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

C  RETVAL=NXP$GETATOMID(,INDIUM_CALIBRATION_PEAK.END\ 

C  1  IDFLAG2,NXP$_ATYPE_SLOT)  !  get  the  R-edge  address 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

RETVAL=NXP$VOLUNTEER(IDFLAG,NXP$_DESC_FLOAT,AREA, 

1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$__ERR_NOERR)  GO  TO  5000 

c  volunteer  the  l-edge  value 

C  RETVAUNXP$VOLUNTEER(IDFLAG1  ,NXP$_DESC_FLOAT,L_EDGE, 

C  1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

C  volunteer  the  r-edge  value 

C  RETVAL=NXP$VOLUNTEER(IDFLAG2,NXP$_DESC_FLOAT,R_EDGE, 

C  1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

GO  TO  95 
ENDIF 


IF  (CODE.EQ.'GA')  THEN  I  if  ifs  a  gallium  peak 
C  get  the  area  address 

RETVAL=NXP$GETATOMID(GALLIUM_CALIBRATION_PEAK.AREA\ 

1  IDFLAG,NXP$_ATYPE_SLOT) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

C  RETVAL=NXP$GETATOMID(,GALLIUM_CALIBRATION_PEAK.START\ 

C  1  IDFLAG1,NXP$_ATYPE_SLOT)  !  get  the  l-edge  address 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

C  RETVAL=NXP$GETATOMID(GALLIUM_CALIBRATION_PEAK.END\ 

C  1  I DF LAG2 , NX P$_ATY P E_S LOT)  I  get  the  R-edge  address 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

C  volunteer  the  area  value 

RETVAL=NXP$VOLUNTEER(IDFLAG,NXP$_DESC_FLOAT,AREA, 

1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

c  volunteer  the  l-edge  value 

C  RETVAL=NXP$VOLUNTEER(!DFLAGl  ,NXP$_DESC_FLOAT,L_EDGE, 

C  1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 
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IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 


C  volunteer  the  r-edge  value 

C  RETVAL=NXP$VOLUNTEER(IDFLAG2,NXP$_DESC_FLOAT,R_EDGE, 

C  1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

C  IF  (RETVAL  .NE.  NXP$_ERR_NOERR)  GO  TO  5000 

GO  TO  95 
ENDIF 

5000  FLAG  =  . FALSE.  I  only  exit  here  if  an  error  occurred 

RETVAL=NXP$GETATOMID(,STATUS.FINDAREA\IDFLAG,NXP$_ATYPE_SLOT) 

RETVAL=NXP$VOLUNTEER(IDFLAG,NXP$_DESC_INT,FLAG, 

1  JIOR(NXP$_VSTRAT_NOFWRD,NXP$M_VSTRAT_SET)) 

9  5  CLOSE(  UNIT=1  )  I  normal  exit  point 
9  6  RETURN 

END 


FUNCTION  IFIND(POINT) 

This  function  returns  the  x  index  value  of  POINT  in  SPECTRA 

COMMON  /BLOCK1/X(20000),  Y(20000),LOADED.NUMBPTS 


I  =  1 

DO  WHILE  (X(I).LT. POINT) 
1  =  1  +  1 

BID  DO 
IFIND=I 
RETURN 
END 


FUNCTION  SLOPE(J) 

This  routine  calculates  the  average  slope  of  a  curve  segment.  It 
currently  uses  10  points  to  average.  The  parameter  J  is  the  index 
number  of  the  leftmost  point  in  the  slope  averaging  ’window’. 


COMMON  /BLOCK1/X(20000),  Y(20000), LOADED, NUMBPTS 

REAL  M,NUM 

M-0.0 

DO  10  1=1  ,  10 

NUM=  Y(J+1)-  Y(J) 

DENOM=  X(J+1)-  X(J) 

IF(DENOM.EQ.O)  DENOM=OLDENOM 
M=  NUM/DENOM  +  M 
OLDENOM=  DENOM 
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J=  J  + 1 
1 0  CONTINUE 

SLOPE-  M/10.0 

RETURN 

0JD 


C . 

Subroutine  PKAREA(STP,  ENP,  AREA) 


This  program,  given  a  set  of  endpoints  and  the  (x,y)  points 
in  between,  will  calculate  the  area  under  the  peak  using 
a  trapezoidal  approximation.  First  calculating  the  area 
from  the  peak  to  thge  zero  baseline  and  then  subtracting 
from  that  the  area  from  the  zero  baseline  to  the  peak  baseline. 

Variable  Declarations 


STP  ->  starting  point  of  the  peak 
ENP  ->  ending  point  of  the  peak 
SPECTRA  ->  all  the  points  in  the  spectra  represented  by 
their  (x,y)  values  in  the  array  form: 

(data  #,  x  or  y)  =  array  representation 

(1.x)  =  (1.1) 

(l.y)  =  (1.2) 

(2.x)  =  (2,1) 

(2,y)  =  (2,2) 


(n,x)  =  (n,1) 

(n,y)  =  (n,2) 

RESAR  ->  area  from  peak  baseline  to  zero  baseline 
SUM  ->  area  under  the  peak  and  above  the  zero  baseline 
AREA  ->  the  area  under  the  peak  and  above  the  peak  baseline 
I  ->  loop  control  variables 

Common  /BLOCK1/X(20000),  Y(20000), LOADED, NUMBPTS 
Real  EP(2),  SP(2) 


Loop  to  find  the  starting  point  in  the  spectra 


I  =  1 

Do  While  (STP  .gt.  X(l)) 

1  =  1  +  1 

End  Do 

1  =  1-1  I  back  up  a  hair 

SP(1)  =  X(l) 

SP(2)  =  Y(l) 
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*Loop  to  calculate  the  area  from  the  peak  to  the  zero  baseline  * 


SUM  -  0.0 

Do  While  (ENP  .ge.  X(l)) 

SUM  =  SUM  +  0.5  *  (X<!+1)  -  X(l))  * 
&  (Y(l)  +  Y(l+1)) 

1-1  +  1 


End  do 

EP(1)  -  X(l) 
EP(2)  =  Y(l) 


Calculate  RESAR  and  subtract  it  from  SUM 


RESAR  .  0.5  *  (EP(1)  -  SP(1))*(EP(2)  +  SP(2)) 
AREA  =  SUM  -  RESAR 
RETURN 
BO 
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APPENDIX  B.  LIST  OF  RULES 


RULE  :  Rule  INVESTIGATE_ALUMINUM_ABS  (#1) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  SPECTRA.TYPE  is  "ABS" 

Then  ALUMINUM_ANALYZED_ABS 
is  confirmed. 

And  CODE  is  set  to  "AL" 

And  Execute 

"FINDAREA"(@ATOMID=SPECTRA.POINTS,ALUMINUM_CALIBRATION_ 

PEAK.WAVENUMBER.CODE.Value;) 

And  (ALUMINUM_CALIBRATlON_PEAK.AREA*4.1 666e+1 3)  is 
assigned  to  ALUMINUM_CALIBRATION_PEAK.CONC 

RULE  :  Rule  INVESTIGATE_ALUMINUM_PHOTO  (#2) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  SPECTRA.TYPE  is  "PHOTO" 

Then  ALUMINUM_ANALYZED_PHOTO 
is  confirmed. 

And  CODE  is  set  to  "AL" 

And  0.0  is  assigned  to  ALUMINUM_CALIBRATION_PEAK.AREA 

RULE  :  Rule  ALUMINUMJJNE_X4  (#3) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  greater 
than  413.5 

And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  less  than 
414.5 

Then  ALUMINUM_LINE 
is  confirmed. 

And  AL_LINE_X4. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  assigned 
to  AL  LINE  X4.WN 


RULE  :  Rule  ALUMINUM_LINE_X3  (#4) 

If 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  greater 
than  403.5 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 

404.5 

Then  ALUMINUM_LINE 
is  confirmed. 

And  AL_LINE_X3. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  AL_LINE_X3.WN 

RULE  :  Rule  ALUMINUM_LINE_X2  (#5) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <[POSSIBLE_ALUMINUM_PEAKSj>.  WAVENUMBER  is  greater 
than  371.5 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 

372.5 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_X2. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigntd 
to  ALJJNE_X2,WN 

RULE  :  Rule  ALUMINUM_LINE_X1  (#6) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  t|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  greater 
than  337.5 

And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  less  than 

338.5 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_X1. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS!>. WAVENUMBER  is  assigned 
to  AL  LINE  XI  .WN 


RULE  :  Rule  ALUMINUMJJNEJEX  (#7) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  greater 
than  453.5 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 
4454.5 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_EX. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  AL_LINE_EX.WN 

RULE  :  Rule  ALUMINUMJJNE  JESTAR  (#8) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  greater 
than  911.5 

And  <|POSSIBLE_ALUMINUM_PEAKS[>. WAVENUMBER  is  less  than 
912 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_ESTAR. EXIST  is  set  to  TRUE 
And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  assigned 
to  AL_LINE_ESTAR.WN 

RULE  :  Rule  ALUMINUM  JJNEJEI  (#9) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  greater 
than  556.5 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 
557 

Then  ALUMINUM_LINE 
is  confirmed. 

And  AL_LINE_EI. EXIST  is  set  to  TRUE 

And  <!POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  assigned 
to  AL  LINE  EI.WN 
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RULE  :  Rule  ALUMINUM_LINE_9  (#10) 

If 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  greater 
than  546 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 
548 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_9. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  AL_LINE_9.WN 

RULE  :  Rule  ALUMINUM_LINE_8  (#11) 

If 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>.  WAVENUMBER  is  greater 
than  543 

And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  less  than 
546 

Then  ALUMINUM_LINE 
is  confirmed. 

And  AL_LINE_8. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  ALJLINE8.WN 

RULE  :  Rule  ALUMINUM_LINE_7  (#12) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  greater 
than  540.7 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 
542 

Then  ALUMINUM_LINE 
is  confirmed. 

And  AL_LINE_7. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  assigned 
to  AL  LINE  7.WN 
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RULE  :  Rule  ALUMINUM_LINE_6  (#13) 

If 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKSj>.WAVENUMBER  is  greater 
than  538 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.  WAVENUMBER  is  less  than 
539.3 

Then  ALUMINUM_LINE 
is  confirmed. 

And  AL_LINE_6. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  AL_LINE_6.WN 

RULE  :  Rule  ALUMINUM_LINE_5  (#14) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>.  WAVENUMBER  is  greater 
than  534.5 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 

535.5 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL__LINE_5. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  AL_LINE_5.WN 

RULE  :  Rule  ALUMINUM_LINE_4P  (#15) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  greater 
than  900 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 

900.5 

Then  ALUMINUM_LINE 
is  confirmed. 

And  AL_LINE_4P. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  AL  LINE  4P.WN 
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RULE  :  Rule  ALUMINUM_L1NE_3P  (#16) 

If 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  greater 
than  891.5 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 
892 

Then  ALUMINUM_LINE 
is  confirmed. 

And  ALJJNEJ3P.  EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  ALJJNEJ3P.WN 

RULE  :  Rule  ALUMINUM_LINE_2P  (#17) 

If 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS(>.  WAVENUMBER  is  greater 
than  867 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 
868 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_2P. EXIST  is  set  to  TRUE 

And  <(POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  assigned 
to  AL_LINE_2P.WN 

RULE  :  Rule  ALUMINUM_LINE_2  (#18) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  greater 
than  471 

And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  less  than 
472.5 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_2. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  AL  LINE  2.WN 
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RULE:  Rule  ALUMINUM  JJNE_1 1  (#19) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  greater 
than  554.5 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 

556.4 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_11. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  AL_LINE_1 1  .WN 

RULE  :  Rule  ALUM1NUM_LINE_10  (#20) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  greater 
than  549 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 
554 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_1 0. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  assigned 
to  AL_LINE_10.WN 

RULE:  Rule  ALUMINUM_LINE_1  (#21) 

I  f 

there  is  evidence  of  ALUMINUM_PRESENT 
And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  greater 
than  442 

And  <|POSSIBLE_ALUMINUM_PEAKS|>.WAVENUMBER  is  less  than 

443.5 

Then  ALUMINUMJJNE 
is  confirmed. 

And  AL_LINE_1.  EXIST  is  set  to  TRUE 

And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  assigned 
to  AL  LINE  1  .WN 


RULE  :  Rule  INVESTIGATE_ALUMINUM  (#22) 

If 

there  is  evidence  of  POSSIBLE_ALUMINUM 
And  <|POSSIBLE_ALUMINUM_PEAKS|>. WAVENUMBER  is  greater 
than  441.7 

And  <jPOSSIBLE_ALUMINUM_PEAKSj>. WAVENUMBER  is  less  than 
443.7 

And  «|POSSIBLE_ALUMINUM_PEAKS|».WAVENUMBER  is  greater 
than  470.9 

And  «|POSSIBLE_ALUMINUM_PEAKS|».WAVENUMBER  is  less 
than  472.9 

And  «<|POSSIBLE_ALUMINUM_PEAKSj»>.WAVENUMBER  is 
greater  than  542.9 

And  «<|POSSIBLE_ALUMINUM_PEAKS|»>.WAVENUMBER  is  less 
than  544.9 

Then  ALUMINUM_PRESENT 
is  confirmed. 

And  «|POSSIBLE_ALUMINUM__PEAKS|». WAVENUMBER  is  assigned 
to  ALUMINUM_CALIBRATION_PEAK.WAVENUMBER 

RULE  :  Rule  START_THINGS_OFF  (#23) 

If 

there  is  evidence  of  PEAKS_LOADED 
And  there  is  evidence  of  PEAKS_EXAMINED 
And  there  is  evidence  of  PEAKS_ANALYZED 
Then  ANALYZED 

is  confirmed. 

RULE  :  Rule  INVESTIGATE_BORON_ABS  (#24) 

If 

there  is  evidence  of  BORON_PRESENT 
And  SPECTRA. TYPE  is  "ABS" 

Then  BORON_AN ALYZE  D_ABS 
is  confirmed. 

And  CODE  is  set  to  "BO" 

And  Execute 

"FINDAREA"(@ATOMID=SPECTRA.POINTS,BORON_CALIBRATION_PEAK. 

WAVENUMBER.CODE.Value;) 

And  (BORON_CALIBRATION_PEAK.AREA*1 .51 51  e+1 3)  is 
assigned  to  BORON_CAUBRATION_PEAK.CONC 
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RULE  :  Rule  INVESTIGATE_BORON_PHOTO  (#25) 

I  f 

there  is  evidence  of  BORON_PRESENT 
And  SPECTRA.TYPE  is  "PHOTO" 

Then  BORON_ANALYZED_PHOTO 
is  confirmed. 

And  CODE  is  set  to  "BO" 

And  0.0  is  assigned  to  BORON_CALIBRATION_PEAK.AREA 

RULE  :  Rule  BORON_LINE_X4  (#26) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  greater  than 
260.0 

And  <|POSSIBLE_BORON_PEAKS(>.WAVENUMBER  is  less  than 
262.0 

Then  BORON _UNE 
is  confirmed. 

And  BO_LINE_X4. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BO_LINE_X4.WN 

RULE  :  Rule  BORON J_INE_X3  (#27) 

I  f 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  greater  than 
249.0 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 
251.0 

Then  BORON_LINE 
is  confirmed. 

And  BO_LINE_X3. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  assigned  to 
BO_LINE_X3.WN 
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RULE  :  Rule  BORON_LINE_X2  (#28) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>.  WAVENUMBER  is  greater  than 
219.0 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  less  than 
221.0 

Then  BORON_LINE 
is  confirmed. 

And  BO_JJNE_X2. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  assigned  to 
BO_LINE_X2.WN 

RULE  :  Rule  BORON_LINE_X1  (#29) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE__BORON_PEAKS|>.WAVENUMBER  is  greater  than 

183.5 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  less  than 

184.5 

Then  BORON_LINE 
is  confirmed. 

And  BO_UNE_X1. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  assigned  to 
BO_LINE_X1.WN 

RULE  :  Rule  BORON_LINE_EX  (#30) 

I  f 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  greater  than 
298.0 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 
300.0 

Then  BORONJJNE 
is  confirmed. 

And  BO_LINE_EX. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BO_LINE_EX.WN 
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RULE  :  Rule  BORONJJN E_EST AR  (#31) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  greater  than 

712.4 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 

713.4 

Then  BORONJJNE 
is  confirmed. 

And  BO_LINE_ESTAR. EXIST  is  set  to  TRUE 
And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BO_LINE_ESTAR.WN 

RULE  :  Rule  BORONJJN E_E  (#32) 

I  f 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  greater  than 
358.0 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 
359.0 

Then  BORONJJNE 
is  confirmed. 

And  BO J.INEJE.  EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  assigned  to 
BO_LINE_E.WN 

RULE  :  Rule  BORON JJNE_9  (#33) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  greater  than 
348.6 

And  <|POSSIBLE_BORON_PEAKSj>. WAVENUMBER  is  less  than 
349.3 

Then  BORONJJNE 
is  confirmed. 

And  BO_LINE_9. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BO_LINE_9.WN 
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RULE  :  Rule  BORON_LINE_8  (#34) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <| POSSI BLE_BORON_P EAKS | > .  WA VE  NU  M B E R  is  greater  than 
344.2 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 
346.7 

Then  BORONJJNE 
is  confirmed. 

And  BO_LINE_8. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  assigned  to 
BO_LINE_8.WN 

RULE  :  Rule  BORON_LINE_7  (#35) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  greater  than 

341.6 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  less  than 

342.6 

Then  BORONJJNE 
is  confirmed. 

And  BOJ_INE_7.EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BO_LINE_7.WN 

RULE  :  Rule  BORON_LINE_6  (#36) 

I  f 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKSj>. WAVENUMBER  is  greater  than 

337.5 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  less  than 

340.6 

Then  BORONJJNE 
is  confirmed. 

And  BO_LINE_6. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  assigned  to 
BO_LINE_6.WN 
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RULE  :  Rule  B0R0N_LINE_5  (#37) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <|POSS!BLE_BORON_PEAKS|>. WAVENUMBER  is  greater  than 
334.0 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 
335.0 

Then  BORONJJNE 
is  confirmed. 

And  BO_LINE_5. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BO_LINE_5.WN 

RULE  :  Rule  BORON  JJNE_4P  (#38) 

I  f 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  greater  than 

700.8 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 

701.8 

Then  BORONJJNE 
is  confirmed. 

And  BO_LINE_4P. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BOJJNE_4P.WN 

RULE  :  Rule  BORON_LINE_4  (#39) 

I  f 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  greater  than 

318.8 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 
322.4 

Then  BORONJJNE 
is  confirmed. 

And  BO_LINE_4. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  assigned  to 
BO  LINE  4.WN 
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RULE  :  Rule  BORON_LINE_3P  (#40) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <iPOSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  greater  than 

692.3 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 

693.4 

Then  BORON_LINE 
is  confirmed. 

And  BO_LINE_3P. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BO_LINE_3P.WN 

RULE  :  Rule  BORON_LINE_3  (#41) 

I  f 

there  is  evidence  of  BORON_P RESENT 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  greater  than 

308.8 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 

309.8 

Then  BORON_LINE 
is  confirmed. 

And  BO_LINE_3. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BOJLINE_3.WN 

RULE  :  Rule  BORON_LINE_2P  (#42) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  greater  than 
668.1 

And  <|POSSlBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 
669.1 

Then  BORON  LINE 
is  confirmed. 

And  BO_LINE_2P. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKSj>.WAVENUMBER  is  assigned  to 
BO_LINE_2P.WN 
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RULE  :  Rule  BORON_LINE_2  (#43) 

I  f 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>.  WAVENUMBER  is  greater  than 

277.7 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 

278.7 

Then  BORONJ.INE 
is  confirmed. 

And  BO_LINE_2. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.  WAVENUMBER  is  assigned  to 
BO_LINE_2.WN 

RULE  :  Rule  BORON_LINE_1 1  (#44) 

I  f 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  greater  than 
356.0 

And  <|POSSIBLE_BORON_PEAKSi>.WAVENUMBER  is  less  than 
357.6 

Then  BORONUNE 
is  confirmed. 

And  BO_LINE_11. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BO_LINE_1 1  .WN 

RULE  :  Rule  BORON_LINE_10  (#45) 

I  f 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  greater  than 
350.3 

And  <|POSSIBLE_BORON_PEAKSj>.WAVENUMBER  is  less  than 
356.0 

Then  BORON_LINE 
is  confirmed. 

And  BO_LINE_10. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  assigned  to 
BO_LINE_10.WN 
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RULE  :  Rule  BORON_LINE_1  (#46) 

If 

there  is  evidence  of  BORON_PRESENT 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  greater  than 

244.5 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 

245.5 

Then  BORON_LINE 
is  confirmed. 

And  BOJJNEJ. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_BORON_PEAKS|>. WAVENUMBER  is  assigned  to 
BOJJNEJ  .WN 

RULE  :  Rule  INVESTIGATE J30R0N  (#47) 

I  f 

there  is  evidence  of  POSSIBLEJBORON 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  greater  than 

244.5 

And  <|POSSIBLE_BORON_PEAKS|>.WAVENUMBER  is  less  than 

245.5 

And  «|POSSIBLE_BORON_PEAKS|».WAVENUMBER  is  greater  than 

277.2 

And  «|POSSIBLE_BORON_PEAKS|».WAVENUMBER  is  less  than 

279.2 

And  «<|POSSIBLE_BORON_PEAKS|»>.WAVENUMBER  is  greater 
than  318.3 

And  «<|POSSIBLE_BORON_PEAKS|»>.WAVENUMBER  is  less  than 

320.3 

Then  BORON_PRESENT 
is  confirmed. 

And  «|POSSIBLE_BORON_PEAKS|».WAVENUMBER  is  assigned  to 
BORON_CALIBRATION_PEAK.  WAVENUMBER 
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RULE  :  Rule  INVESTIGATE_CARBON_ABS  (#48) 

If 

there  is  evidence  of  CARBON_PRESENT 
And  SPECTRA.TYPE  is  HABSH 
Then  CARBONANALYZED_ABS 
is  confirmed. 

And  CODE  is  set  to  "CA" 

And  Fxecute 

"FINDAREA"(@ATOMID=SPECTRA.POINTS,CARBON_CALIBRATION_PEAK. 
WAVENUMBER, CODE.Value;) 

And  (CARBON_CALIBRATION_PEAK.AREA*6.700e+16)  is 
assigned  to  CARBON_CALIBRATION_PEAK.CONC 

RULE  :  Rule  INVESTIGATE_CARBON_PHOTO  (#49) 

I  f 

there  is  evidence  of  CARBON_PRESENT 
And  SPECTRA.TYPE  is  "PHOTO” 

Then  CARBON__ANALYZED_PHOTO 
is  confirmed. 

And  CODE  is  set  to  "CA" 

And  0.0  is  assigned  to  CARBON_CALIBRATION_PEAK.AREA 

RULE  :  Rule  CARBON_LINE_1  (#50) 

I  f 

there  is  evidence  of  CARBON_PRESENT 

And  <|POSSIBLE_CARBON_PEAKS|>.WAVENUMBER  is  greater  than 

599.5 

And  <|POSSIBLE_CARBON_PEAKSj>.WAVENUMBER  is  less  than 

600.5 

Then  CARBONJJNE 
is  confirmed. 

And  CA_LINE_1. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_CARBON_PEAKS|>.WAVENUMBER  is  assigned  to 
CA  LINE  1.WN 


RULE  :  Rule  INVESTIGATE_CARBON  (#51) 

If 

there  is  evidence  of  POSSIBLE_CARBON 

And  <|POSSIBLE_CARBON_PEAKS|>.WAVENUMBER  is  greater  than 

599.5 

And  <|POSSIBLE_CARBON_PEAKS|>.WAVENUMBER  is  less  than 

600.5 

Then  CARBON_PRESENT 
is  confirmed. 

And  <|POSSIBLE_CARBON_PEAKS|>.WAVENUMBER  is  assigned  to 
CARBON_CALIBRATION_PEAK.WAVENUMBER 

RULE  :  Rule  ALJJNES  (#52) 

I  f 

there  is  evidence  of  ANALYZED 
And  there  is  evidence  of  ALUMINUM_PRESENT 
Then  FIND_AL_LINES 
is  confirmed. 

And  ALUMINUM_LINE  is  assigned  to  ALUMINUMJJNE 

RULE  :  Rule  BO_LINES  (#53) 

If 

there  is  evidence  of  ANALYZED 
And  there  is  evidence  of  BORON_PRESENT 
Then  FIND_BO_LINES 
is  confirmed. 

And  BORON_LINE  is  assigned  to  BORONJJNE 

RULE  :  Rule  CAJJNES  (#54) 

If 

there  is  evidence  of  ANALYZED 
And  there  is  evidence  of  CARBON_PRESENT 
Then  FIND_CA_LINES 
is  confirmed. 

And  CARBONJJNE  is  assigned  to  CARBON_LINE 
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RULE  :  Rule  GAJJNES  (#55) 

I  f 

there  is  evidence  of  ANALYZED 
And  there  is  evidence  of  GALLIUM_PRESENT 
Then  FIND_GA_LINES 
is  confirmed. 

And  GALLIUMJJNE  is  assigned  to  GALLIUMJJNE 

RULE  :  Rule  INJJNES  (#56) 

I  f 

there  is  evidence  of  ANALYZED 
And  there  is  evidence  of  INDIUM_PRESENT 
Then  FINDJNJJNES 
is  confirmed. 

And  INDIUMJJNE  is  assigned  to  INDIUMJJNE 

RULE  :  Rule  OX_LINES  (#57) 

I  f 

there  is  evidence  of  ANALYZED 
And  there  is  evidence  of  OXYGEN_PRESENT 
Then  FIND_OX_LINES 
is  confirmed. 

And  OXYGENJJNE  is  assigned  to  OXYGEN_LINE 

RULE  :  Rule  INVESTIGATE_GALLIUM_ABS  (#58) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  SPECTRA.TYPE  is  "ABS" 

Then  GALLIUM_ANALYZED_ABS 
is  confirmed. 

And  CODE  is  set  to  "GA" 

And  Execute 

"FINDAREA"(@ATOMID=SPECTRA.POINTS,GALLIUM_CALIBRATION_PEAK 

.WAVENUMBER.CODE.Value;) 

And  (GALLIUM_CALIBRATION_PEAK.AREA*2.780e+14) 
is  assigned  to  GALLIUM_CALIBRATION_PEAK.CONC 
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RULE  :  Rule  INVESTIGATE_GALLIUM_PHOTO  (#59) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  SPECTRA.TYPE  is  "PHOTO" 

Then  GALLIUM_ANALYZED_PHOTO 
is  confirmed. 

And  CODE  is  set  to  "GA" 

And  0.0  is  assigned  to  GALLIUM_CALIBRATION_PEAK.AREA 

RULE  :  Rule  GALLIUM_LINE_X4  (#60) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  greater 
than  422 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
424 

Then  GALLIUMJJNE 
is  confirmed. 

And  GAJJNE_X4. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
GA_LINE_X4.WN 

RULE  :  Rule  GALLIUM_LINE_X2  (#61) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  380 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
382 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_X2. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_X2.WN 
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RULE  :  Rule  GALLIUM_LINE_X1  (#62) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  346 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  less  than 
348 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_X1  .EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_X1.WN 

RULE  :  Rule  GALLIUM_LINE_EX  (#63) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS[>. WAVENUMBER  is  greater 
than  459 

And  <|POSSIBLE_GALL!UM_PEAKS|>. WAVENUMBER  is  less  than 
461 

Then  GALLIUM_LINE 
is  confirmed. 

And  GA_LINE_EX. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_EX.WN 

RULE  :  Rule  GALLIUM_LINE_ESTAR  (#64) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  941 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
942 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_ESTAR. EXIST  is  set  to  TRUE 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
GA  LINE  ESTAR.WN 
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RULE  :  Rule  GALLIUM_LINE_EI  (#65) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  586 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  less  than 

587.2 

Then  GALLIUMJJNE 
is  confirmed. 

And  GAJJNE_EI. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_EI.WN 

RULE  :  Rule  GALLIUM_LINE_9  (#66) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  576 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 

578.2 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_9. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_9.WN 

RULE  :  Rule  GALLIUM  JJNEJJ  (#67) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  greater 
than  573 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
576 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_8. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GAJINE_8.WN 
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RULE  :  Rule  GALL!UM_LINE_7  (#68) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  570.5 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 

571.5 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_7. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GAJJNE_7.WN 

RULE  :  Rule  GALLIUM_LINE_6  (#69) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  567 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 

569.5 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_6. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALL!UM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_6.WN 

RULE  :  Rule  GALLIUM  JJNEJ5P  (#70) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  greater 
than  933.5 

And  <|POSSIBLE_GALLIUM_PEAKSj>.WAVENUMBER  is  less  than 

934.5 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_5P. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_5F.WN 
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RULE  :  Rule  GALLIUM_LINE_5  (#71) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  greater 
than  563.7 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
564.7 

Then  GALLIUM_LINE 
is  confirmed. 

And  GAJJNE_5. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_5.WN 

RULE  :  Rule  GALLIUM_LINE_4P  (#72) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  929.5 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  less  than 
930.5 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_4P. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_4P.WN 

RULE  :  Rule  GALLIUM_LINE_4  (#73) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  < | POSS I BLE_G ALU UM_PEAKS|>.WAVENUMBER  is  greater 
than  547.7 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
553 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_4. EXIST  is  set  to  TRUE 

And  <jPOSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA_LINE_4.WN 
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RULE  :  Rule  GALLIUM_LINE_3P  (#74) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  921 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
922 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_3P. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
GA_LINE_3P.WN 

RULE  :  Rule  GALLIUM_LINE_3  (#75) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  541 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  less  than 
542 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_L!NE_3. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
GA_LINE_3.WN 

RULE  :  Rule  GALLIUM JJNE_2P  (#76) 

I  f 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  greater 
than  896.7 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
897.7 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_2P. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKSj>. WAVENUMBER  is  assigned  to 
GA  LINE  2P.WN 
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RULE  :  Rule  GALLIUM _LINE_1 1  (#77) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  < jPOSSIBLE_GALLIUM_PEAKS| >.WAVENUMBE R  is  greater 
than  584 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  less  than 
586.3 

Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_11. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLiUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GAJJNEJ 1  .WN 

RULE  :  Rule  GALLIUM_LINE_1 0  (#78) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  579 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
584 

Then  GALLIUM_LINE 
is  confirmed. 

And  GAJJNE_1 0. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
GAJJNEJ 0.WN 

RULE  :  Rule  GALLIUM_LINE_1  (#79) 

If 

there  is  evidence  of  GALLIUM_PRESENT 
And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  greater 
than  469.25 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  less  than 
470.25 
Then  GALLIUMJJNE 
is  confirmed. 

And  GA_LINE_1. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_GALLIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
GA  LINE  1  .WN 
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RULE  :  Rule  INVESTIGATE_GALLIUM  (#80) 

I  f 

there  is  evidence  of  POSSIBLE_GALLIUM 
And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  greater 
than  468.8 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  less  than 
470.8 

And  «|POSSIBLE_GALLIUM_PEAKS|».WAVENUMBER  is  greater 
than  540.5 

And  «|POSSIBLE_GALLIUM_PEAKS|».WAVENUMBER  is  less  than 
542.5 

And  «<|POSSIBLE_GALLIUM_PEAKS|»>.WAVENUMBER  is  greater 
than  547.2 

And  «<|POSSIBLE_GALLIUM_PEAKS|»>.WAVENUMBER  is  less 
than  549.2 

Then  GALLIUM_PRESENT 
is  confirmed. 

And  <|POSSIBLE_GALLIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
GALLIUM_CALIBRATION_PEAK.WAVENUMBER 

RULE  :  Rule  INVESTIGATE_INDIUM_ABS  (#81) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 
And  SPECTRA.TYPE  is  "ABS" 

Then  INDIUM_ANALYZED_ABS 
is  confirmed. 

And  CODE  is  set  to  "IN" 

And  Execute 

"FINDAREA"«3>ATOMID=SPECTRA.POINTS,INDIUM_CALIBRATION_PEAK. 
WAVENUMBER, CODE.Value;) 

And  (INDIUM_CALIBRATION_PEAK.AREA*2.0000000000e+1 5)  is 
assigned  to  INDIUM_CALIBRATION_PEAK.CONC 

RULE  :  Rule  INVESTIGATE_INDIUM_PHOTO  (#82) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 
And  SPECTRA.TYPE  is  "PHOTO" 

Then  INDIUM_ANALYZED_PHOTO 
is  confirmed. 

And  CODE  is  set  to  "IN” 

And  0.0  is  assigned  to  INDIUM_CALIBRATlON_PEAK.AREA 
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RULE  :  Rule  INDIUM_LINE_X4  (#83) 

If 

there  is  evidence  of  IND1UM_PRESENT 

And  <|POSSIBLE_INDIUM  PEAKS|>.WAVENUMBER  is  greater  tf  an 
872 

And  <|POSSIBLEJNDIUM_PEAKSj>. WAVENUMBER  is  less  than  874 
Then  INDIUMJJNE 
is  confirmed. 

And  IN_LINE_X4.EXIST  is  set  to  TRUE 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
INJ.INE_X4.WN 

RULE  :  Rule  INDIUM  JJNE_X3  (#84) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  greater  than 
862 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than  864 
Then  INDIUMJJNE 
is  confirmed. 

And  IV  'JNE_X3.EXIST  is  set  to  TRUE 

And  <iPOSSIBLEJNDIUM_PEAKSj>. WAVENUMBER  is  assigned  to 
IN_LINE_X3.WN 

RULE  :  Rule  INDIUM_LINE_X2  (#85) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  greater  than 
830 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  less  than  832 
Then  INDIUMJJNE 
is  confirmed. 

And  IN  JJNE_X2. EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN  LINE  X2.WN 
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RULE  :  Rule  INDIUM_LINE_X1  (#86) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  greater  than 
798 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than  780 
Then  INDIUMJJNE 
is  confirmed. 

And  IN_LINE_X1. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
IN_LINE_X1  .WN 

RULE  :  Rule  INDIUM_LINE_EX  (#87) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  greater  than 
908 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  less  than  910 
Then  INDIUM_LINE 
is  confirmed. 

And  IN_LJNE_EX.  EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN_LINE_EX.WN 

RULE  :  Rule  INDIUM__LINE_ESTAR  (#88) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>.  WAVENUMBER  is  greater  than 

1609.2 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than 

1610.2 
Then  INDIUMJJNE 

is  confirmed. 

And  IN_LINE_ESTAR. EXIST  is  set  to  TRUE 
And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN  LINE  ESTAR.WN 
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RULE  :  Rule  INDIUM_LINE_EI  (#89) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  greater  than 

1254.3 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  less  than 

1255.3 
Then  INDIUM_LINE 

is  confirmed. 

And  IN_LINE_EI. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_INDIUM_PEAKSj>. WAVENUMBER  is  assigned  to 
IN_LINE_EI.WN 

RULE  :  Rule  INDIUM_LINE_9  (#90) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  greater  than 

1245.3 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than 

1246.3 
Then  INDIUM_LINE 

is  confirmed. 

And  IN_LINE_9. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN_LINE_9.WN 

RULE:  Rule  INDIUM_LINE_8  (#91) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  greater  than 
1241 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  less  than 
1244 

Then  INDIUM_LINE 
is  confirmed. 

And  IN_LINE_8. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN  LINE  8.WN 
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RULE  :  Rule  INDIUM_L!NE_7  (#92) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLEJNDIUM_PEAKS|>.WAVENUMRER  is  greater  than 

1238.5 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than 

1239.5 
Then  INDIUMJJNE 

is  confirmed. 

And  IN_LINE_7. EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
IN_LINE_7.WN 

RULE  :  Rule  INDIUM_LINE_6  (#93) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  greater  than 

1236 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than 

1237 

Then  INDIUMJJNE 
is  confirmed. 

And  IN_LINE_6. EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
INJJNE_6.WN 

RULE  :  Rule  INDIUM_LINE_5  (#94) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLEJNDlUM_PEAKS|>. WAVENUMBER  is  greater  than 

1231.7 

And  <|POSSIBLEJNDIUM_PEAKSj>. WAVENUMBER  is  less  than 

1232.7 
Then  INDIUMJJNE 

is  confirmed. 

And  INJJNEJ5. EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKSj>. WAVENUMBER  is  assigned  to 
IN  LINE  5.WN 


105 


RULE  :  Rule  INDIUM_LINE_4P  (#95) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  greater  than 

1598 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than 

1599 

Then  INDIUMJJNE 
is  confirmed. 

And  IN_LINE_4P. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN_LINE_4P.WN 

RULE  :  Rule  INDIUM_LINE_4A  (#96) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  greater  than 

1217.5 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  less  than 

1218.5 
Then  INDIUMJJNE 

is  confirmed. 

And  IN JJNE_4A.  EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN_LINE_4A.WN 

RULE  :  Rule  INDIUM JJNE_4  (#97) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  greater  than 

1215.3 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  less  than 

1217.3 
Then  INDIUMJJNE 

is  confirmed. 

And  IN_LINE_4. EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN  LINE  4.WN 
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RULE  :  Rule  INDIUM_LINE_3P  (#98) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  greater  than 

1589.3 

And  <|POSSIBLE_INDIUM_PEAKSj>.WAVENUMBER  is  less  than 

1590.3 
Then  INDIUMJJNE 

is  confirmed. 

And  IN_LINE_3P. EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
IN_LINE_3P.WN 

RULE  :  Rule  INDIUM_LINE_3  (#99) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  greater  than 

1207.5 

And  <|POSSIBLEJNDIUM_PEAKS|>.WAVENUMBER  is  less  than 

1208.5 
Then  INDIUMJJNE 

is  confirmed. 

And  IN JJNE_3. EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN_LINE_3.WN 

RULE:  Rule  INDIUM_LINE_2P  (#100) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLEJNDIUM_PEAKS|>.WAVENUMBER  is  greater  than 
1565 

And  <|POSSIBLEJNDIUM_PEAKS|>.WAVENUMBER  is  less  than 
1567 

Then  INDIUM_LINE 
is  confirmed. 

And  IN_LINE_2P. EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN  LINE  2P.WN 
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RULE  :  Rule  INDIUM_LINE_2  (#101) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLEJNDIUM_PEAKS|>. WAVENUMBER  is  greater  than 

1175.3 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than 

1176.3 
Then  INDIUMJJNE 

is  confirmed. 

And  INJJNE_2.EXIST  is  set  to  TRUE 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  assigned  to 
IN_LINE_2.WN 

RULE  :  Rule  INDIUM_LINE_1 1  (#102) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  greater  than 
1253 

And  <|POSSIBLEJNDIUM_PEAKS|>.WAVENUMBER  is  less  than 
1255 

Then  INDIUM_LINE 
is  confirmed. 

And  IN_LINE_11. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
IN  _LINE_1 1  .WN 

RULE  :  Rule  INDIUM_LINE_10  (#103) 

If 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  greater  than 

1247.4 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than 

1251.4 
Then  INDIUMJJNE 

is  confirmed. 

And  IN JJNE_10. EXIST  is  set  to  TRUE 

And  <|POSSIBLEJNDIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
IN  LINE  10.WN 
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RULE  :  Rule  INDIUM  JJNEJ  (#104) 

I  f 

there  is  evidence  of  INDIUM_PRESENT 

And  <|POSSIBLE_INDIUM_PEAKS|>.  WAVENUMBER  is  greater  than 

1145 

And  <!POSSIBLE_INDIUM_PEAKS|>. WAVENUMBER  is  less  than 

1146 

Then  INDIUMJJNE 
is  confirmed. 

And  IN_LINE_1. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  assigned  to 
IN_LINE_1  .WN 

RULE  :  Rule  INVESTIGATE JNDIUM  (#105) 

I  f 

there  is  evidence  of  POSSIBLEJNDIUM 

And  <|POSSIBLEJNDIUM_PEAKS|>.WAVENUMBER  is  greater  than 

1144.4 

And  <|POSSIBLE_INDIUM_PEAKS|>.WAVENUMBER  is  less  than 

1146.4 

And  «|POSSIBLE_INDIUM_PEAKS|».WAVENUMBER  is  greater 
than  1174.8 

And  «|POSSIBLE_INDIUM_PEAKS|».WAVENUMBER  is  less  than 
1176.8 

And  «<|POSSIBLEJNDIUM_PEAKS|»>.WAVENUMBER  is  greater 
than  1215.3 

And  <«|POSSIBLE_INDIUM_PEAKS|>».WAVENUMBER  is  less  than 
1217.3 

Then  INDIUM_PRESENT 
is  confirmed. 

And  «|POSSIBLE_INDIUM_PEAKS|».WAVENUMBER  is  assigned  to 
INDIUM_CALIBRATION_PEAK.WAVENUMBER 
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RULE  :  Rule  INVESTIGATE_OXYGEN_ABS  (#106) 

If 

there  is  evidence  of  OXYGEN_PRESENT 
And  SPECTRA.TYPE  is  "ABS" 

Then  OXYG  EN_ANALYZ  E  D_ABS 
is  confirmed. 

And  CODE  is  set  to  "OX" 

And  Execute 

"FINDAREA"(@ATOMID=SPECTRA.POINTS,OXYGEN_CALIBRATION_PEAK. 
WAVENUMBER,  CODE.Vaiue;) 

And  (OXYGEN_CALIBRATION_PEAK.AREA*3.090e+16)  is  assigned 
to  OXYGEN_CALIBRATION_PEAK.CONC 

RULE  :  Rule  INVESTIGATE_OXYGEN_PHOTO  (#107) 

I  f 

there  is  evidence  of  OXYGEN_PRESENT 
And  SPECTRA.TYPE  is  "PHOTO" 

Then  OXYGEN_ANALYZEDJPHOTO 
is  confirmed. 

And  CODE  is  set  to  "OX" 

And  0.0  is  assigned  to  OXYGEN_CALIBRATION_PEAK.AREA 

RULE:  Rule  OXYGEN_LINE_1  (#108) 

If 

there  is  evidence  of  OXYGEN_PRESENT 

And  <|POSSIBLE_OXYGEN_PEAKS|>.WAVENUMBER  is  greater  than 
1135.0 

And  <|POSSIBLE_OXYGEN_PEAKS|>.WAVENUMBER  is  less  than 
1137.0 
Then  OXYGEN  LINE 
is  confirmed. 

And  OX_LINE_1. EXIST  is  set  to  TRUE 

And  <|POSSIBLE_OXYGEN_PEAKS|>. WAVENUMBER  is  assigned  to 
OX_LINE_1.WN 
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RULE  :  Rule  INVESTIGATE_OXYGEN  (#109) 

I  f 

there  is  evidence  of  POSSIBLE_OXYGEN 

And  <|POSSIBLE_OXYGEN_PEAKS|>.WAVENUMBER  is  greater  than 
1135.0 

And  <|POSSIBLE_OXYGEN_PEAKS|>.WAVENUMBER  is  less  than 
1137.0 

Then  OXYGEN_PRESENT 
is  confirmed. 

And  <|POSSIBLE_OXYGEN_PEAKS|>.WAVENUMBER  is  assigned  to 
OXYGEN_CALIBRATION_PEAK.WAVENUMBER 

RULE  :  Rule  EXAMINE_THE_POSSIBLE_PEAKS  (#110) 

I  f 

GALLIUM_ANALYZED_ABS  is  not  NOTKNOWN 
And  GALLIUM_ANALYZED_PHOTO  is  not  NOTKNOWN 
And  BORON_ANALYZED_ABS  is  not  NOTKNOWN 
And  BORON_ANALYZED_PHOTO  is  not  NOTKNOWN 
And  INDIUM_ANALYZED_ABS  is  not  NOTKNOWN 
And  INDIUM_ANALYZED_PHOTO  is  not  NOTKNOWN 
And  ALUMINUM_ANALYZED_ABS  is  not  NOTKNOWN 
And  ALUMINUM_ANALYZED_PHOTO  is  not  NOTKNOWN 
And  CARBON_ANALYZED_ABS  is  not  NOTKNOWN 
And  CARBON_ANALYZED_PHOTO  is  not  NOTKNOWN 
And  OXYG EN_ANALYZE D_ABS  is  not  NOTKNOWN 
And  OXYGEN_ANALYZED_PHOTO  is  not  NOTKNOWN 
Then  PEAKS_ANALYZED 
is  confirmed. 

RULE  :  Rule  Force_the_impurity_analysis  (#111) 

I  f 

POSSIBLE_ALUMINUM  is  not  NOTKNOWN 
And  POSSIBLEJ3ALLIUM  is  not  NOTKNOWN 
And  POSSIBLEJNDIUM  is  not  NOTKNOWN 
And  POSSIBLE_BORON  is  not  NOTKNOWN 
And  POSSIBLE_CARBON  is  not  NOTKNOWN 
And  POSSIBLE_OXYGEN  is  not  NOTKNOWN 
And  POSSIBLE_UNKNOWN_PEAKS  is  not  NOTKNOWN 
Then  PEAKS_EXAMINED 
is  confirmed. 


RULE  :  Rule  Load_the_peak_table  (#112) 

I  f 

there  is  evidence  of  SPECTRA_ID_ASKED 
And  SPECTRA. TYPE  is  "ABS", "PHOTO" 

Then  PEAKS_LOADED 
is  confirmed. 

And  Retrieve  from  "@V(SPECTRA.PEAKS)" 
@TYPE«NXPDB;@RLL»ADD;@UNKNOWN=TRUE;@NAME="!NAME!"; 
@CREATE«|POSSIBLE_PEAKS|;@PROPS*AMPLITUDE,WAVENUMBER> 
HALF_WIDTH,AREA;@FIELDS»"AMPLITUDE","WAVENUMBER',> 
"HALF_WIDTH","AREA"; 

RULE  :  Rule  PICK_POSSIBLE_ALUMINUM_PEAKS  (#113) 

If 

<|POSSIBLE_PEAKS|>.WAVENUMBER  is  greater  than  337 
And  <|POSSIBLE__PEAKS|>.WAVENUMBER  is  less  than  912 
Then  POSSIBLE_ALUMlNUM 
is  confirmed. 

And  Create  Object  <|POSSIBLE_PEAKS|> 
|POSSIBLE_ALUMINUM_PEAKS| 

RULE  :  Rule  PICK_OUT_POSSIBLE_BORON_PEAKS  (#114) 

If 

<|POSSIBLE_PEAKS|>.WAVENUMBER  is  greater  than  183 
And  <|POSSIBLE_PEAKS|>.WAVENUMBER  is  less  than  713.5 
Then  POSSIBLE_BORON 
is  confirmed. 

And  Create  Object  <|POSSIBLE_PEAKS|> 
|POSS!BLE_BORON_PEAKS| 

RULE  :  Rule  PICK_POSSIBLE_CARBON_PEAK  (#115) 

If 

<|POSSIBLE_PEAKS|>.WAVENUMBER  is  greater  than  599.0 
And  <|POSSIBLE_PEAKS|>.WAVENUMBER  is  less  than  601 .0 
Then  POSSIBLE_CARBON 
is  confirmed. 

And  Create  Object  <|POSSIBLE_PEAKS|> 
|POSSIBLE_CARBON_PEAKS| 
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RULE  :  Rule  P!CK_POSSIBLE_GALLIUM_PEAKS  (#116) 

I  f 

<|POSSIBLE_PEAKS|>.WAVENUMBER  is  greater  than  346 
And  <|POSSIBLE_PEAKS|>.WAVENUMBER  is  less  than  942 
Then  POSSIBLE_GALLIUM 
is  confirmed. 

And  Create  Object  <|POSSIBLE_PEAKS|> 

|  POSS I BL  E_G  ALLI U  M_P  E  AKS  | 

RULE  :  Rule  PICK_OUT_POSSIBLE_INDIUM_PEAKS  (#117) 

I  f 

<|POSSIBLE_PEAKS|>.WAVENUMBER  is  greater  than  798 
And  <|POSSIBLE_PEAKS|>. WAVENUMBER  is  less  than  1610 
Then  POSSIBLEJNDIUM 
is  confirmed. 

And  Create  Object  <|POSSIBLE_PEAKS|> 

|  POSSIBLE JNDIUM_PEAKS| 

RULE  :  Rule  PICK_POSSIBLE_OXYGEN_PEAKS  (#118) 

If 

<|POSSIBLE_PEAKS|>. WAVENUMBER  is  greater  than  1135.0 
And  <|POSSIBLE_PEAKS|>.WAVENUMBER  is  less  than  1137.0 
Then  POSSIBLE_OXYGEN 
is  confirmed. 

And  Create  Object  <|POSSIBLE_PEAKS|> 
|POSSIBLE_OXYGEN_PEAKS| 


RULE  :  Rule  UNIDENTIFIED_PEAKS  (#119) 

If 

<|POSSIBLE_PEAKS|>  is  not  a  member  of 
<|POSSIBLE_ALUMINUM_PEAKS|> 

And  <|POSSIBLE_PEAKS)>  is  not  a  member  of 
<|POSSIBLE_BORON_PEAKS|> 

And  <|POSSIBLE_PEAKS|>  is  not  a  member  of 
<)POSSIBLE_INDIUM_PEAKS|> 

And  <jPOSSIBLE_PEAKS|>  is  not  a  member  of 
<|POSSIBLE_GALLIUM_PEAKS|> 

And  <|POSSIBLE_PEAKS|>  is  not  a  member  of 
<  |  POSSI  BLE_CARBON_P  EAKS  |  > 

And  <|POSSIBLE_PEAKS|>  is  not  a  member  of 
<|POSSIBLE_OXYGEN_PEAKS|> 

Then  POSSIBLE_UNKNOWN_PEAKS 
is  confirmed. 

And  Create  Object  <|POSSIBLE_PEAKS|>  |UNIDENTIFIED_PEAKS| 

RULE  :  Rule  PLOT_THE_SPECTRA  (#120) 

If 

REPORT_ALUMINUM  is  not  NOTKNOWN 
And  REPORT_BORON  is  not  NOTKNOWN 
And  REPORT_GALLIUM  is  not  NOTKNOWN 
And  REPORTJNDIUM  is  not  NOTKNOWN 
And  REPORTOXYGEN  is  not  NOTKNOWN 
And  REPORT_CARBON  is  not  NOTKNOWN 
Then  reports_done 
is  confirmed. 

And  Execute  "rgraph"(@ATOMID=SPECTRA.NAME;) 

RULE  :  Rule  Get_the_spectra_ID  (#121) 

If 

SPECTRA.NAME  is  assigned  to  SPECTRA.NAME 
Then  SPECTRA_ID_ASKED 
is  confirmed. 

And  Execute  "FILLSPECTRA"(@ATOMID=SPECTRA.NAME;) 

And  Execute  "PEAKPICK"() 

And  STRCAT(SPECTRA. NAME,". MAS")  is  assigned  to 
SPECTRA.POINTS 

And  STRCAT(SPECTRA.NAME,".NXP”)  is  assigned  to 
SPECTRA.  PEAKS 
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